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20.  ABSTRACT  (Continued). 

elastic-plastic,  stress-strain  relationships  are  reviewed  with  respect  to  their 
capabilities  for  modeling  the  stress-strain  behavior  of  compacted  clays  .3  It  is 
shown  that  if  the  Cam  Clay  Model  is  revised  by  introducing  a cohesion  Intercept 
and  a nonlinear  e - log  a ' curve,  it  can  be  used  to  model  the  stress-strain 
behavior  of  compacted  clay  quite  accurately.  This  model  provides  an  incre- 
mental stress-strain  relationship  which  can  be  used  for  analyses  of  drained, 
undrained,  and  partly  drained  conditions . 

Available  data  are  reviewed  to  show  how  the  permeability  of  compacted 
clays  varies  with  degree  of  saturation  and  void  ratio,  and  an  empirical  rela- 
tionship is  developed  which  expresses  the  dependence  of  permeability  on  these 
factors  for  a wide  variety  of  compacted  clays.  A concept  of  homogenized  pore 
fluid  is  proposed  to  simulate  the  air-water  mixture  in  the  pores  of  partly  sat- 
urated clays  so  that  the  three-phase  partly  saturated  clay  can  be  treated 
phenomenologically  as  a two-phase  material.  Compressibility  of  the  homogenized 
pore  fluid  is  then  derived  using  Boyle's  Law  and  Henry's  Law  and  taking  account 
of  the  effect  of  surface  tension.  Biot's  theory  of  consolidation  is  reviewed, 
and  a finite  element  formulation  for  the  consolidation  of  saturated  soils  is 
derived.  The  computer  program  based  on  this  formulation  is  verified  by  com- 
paring results  derived  from  the  program  with  known  solutions  for  several 
problems . 

The  theory  of  consolidation  for  partly  saturated  soil  is  derived  and  a 
computer  program  is  developed.  To  check  the  utility  and  effectiveness  of  the 
theory  and  the  computer  program  developed,  two  examples  of  one-dimensional 
consolidation  with  compressible  pore  fluid  and  varying  permeability  are 
analyzed.  Finite  element  analyses  performed  using  this  program  are  shown  to 
agree  well  with  closed-form  solutions  for  the  consolidation  of  both  saturated 
and  unsaturated  soils.  A finite  element  analysis  is  performed  to  study  the 
.stresses  and  movements  in  New  Melones  Dam,  which  is  now  under  construction. 
^j^The  results  of  this  analysis  indicate  that  the  behavior  of  the  dam  during  con- 
solidation is  closely  related  to  the  stiffness  and  degree  of  saturation  of  the 
core.  The  results  of  this  study  indicate  that  the  finite  element  procedure 
developed  is  a potentially  useful  tool.  It  seems  likely  that  its  greatest 
value  in  application  will  be  in  connection  with  interpretation  of  instrumenta- 
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FOREWORD 

The  work  described  in  this  report  was  performed  under  Contract  No. 
DACW39-74-C-0027,  "Behavior  of  Zoned  Embankments  on  Soft  Foundations," 
between  the  U.  S.  Army  Engineer  Waterways  Experiment  Station  and  the 
University  of  California.  This  is  the  eighth  report  on  investigations  per- 
formed under  this  contract.  The  first  report,  "Finite  Element  Analyses  of 
Stresses  and  Movements  in  Embankments  During  Construction,"  by  F.  H.  Kulhawy, 

J.  M.  Duncan  and  H.  Bolton  Seed,  was  published  in  November,  1969.  The  second 
report,  "Three-Dimensional  Finite  Element  Analyses  of  Dams,"  by  Guy  Lefebvre 
and  J.  M.  Duncan,  was  published  in  May,  1971.  The  third  report,  "Effect  of 
Reservoir  Filling  on  Stresses  and  Movements  in  Earth  and  Rockfill  Dams,"  by 
E.  S.  Nobari  and  J.  M.  Duncan,  was  published  in  January,  1972.  The  fourth 
report,  "Hydraulic  Fracturing  in  Zoned  Earth  and  Rockfill  Dams,"  by  E.  S. 
Nobari,  K.  L.  Lee  and  J.  M.  Duncan,  was  published  in  January,  1973.  The 
fifth  report,  "Finite  Element  Analyses  of  Transverse  Cracking  in  Low  Embank- 
ment Dams,"  by  Guy  Lefebvre  and  J.  M.  Duncan,  was  published  in  October,  1974. 
The  sixth  report,  "Finite  Element  Analyses  of  Stresses  and  Movements  in 
Birch  Dam,"  by  Antonio  Soriano,  J.  M.  Duncan,  Kai  Wong  and  Jean-Michel  Simon, 
was  published  in  April,  1976.  The  seventh  report,  "The  Role  of  Fill  Strength 
in  the  Stability  of  Embankments  on  Soft  Clay  Foundations,"  by  Suphon  Chirapuntu 
and  J.  M.  Duncan,  was  published  in  June,  1976.  The  research  was  sponsored  by 
The  Office,  Chief  of  Engineers,  under  the  Civil  Works  Program  CWIS  No.  31173, 
Special  Studies  for  Civil  Works  Soils  Problems,  Task  2. 

The  general  objective  of  this  research,  which  was  begun  in  June,  1968, 
is  to  develop  methods  for  analysis  of  stresses  and  movements  in  embankments. 
Work  on  this  project  is  conducted  under  the  supervision  of  J.  M.  Duncan, 
Professor  of  Civil  Engineering.  The  project  is  administered  by  the  Office  of 
Research  Services  of  the  College  of  Engineering.  The  phase  of  the  investiga- 
tion described  in  this  report  was  conducted,  and  the  report  was  prepared  by 
C.  S.  Chang  and  J.  M.  Duncan. 

The  contract  was  monitored  by  Mr.  C.  L.  McAnear,  Chief,  Soil  Mechanics 
Division,  under  the  general  supervision  of  Mr.  J.  P.  Sale,  Chief,  Soils  and 

1 Pavements  Laboratory.  Contracting  Officer  was  Col.  John  L.  Cannon,  Director 

of  the  U.  S.  Army  Waterways  Experiment  Station. 
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INTRODUCTION 


INTRODUCTION 

The  theory  of  consolidation  was  first  developed  by  Terzaghi  (1924) . 
In  one-dimensional  cases,  such  as  those  considered  by  Terzaghi,  the  con- 
solidation process  is  governed  by  the  diffusion  equation.  In  two-  or 
three-dimensional  cases  it  is  often  assumed  that  the  consolidation  pro- 
cess is  governed  by  a two-  or  three-dimensional  diffusion  equation,  as, 
for  example,  in  the  work  done  by  Rendulic  (1936) . 

For  problems  with  simple  geometries  and  simple  boundary  condi- 
tions it  is  often  possible  to  find  an  analytic  solution  to  the  diffusion 
equation.  Solutions  of  this  type  have  been  developed  by  Carrilo  (1942) , 
Carslaw  and  Jaegar  (1959) , and  Scott  (1963) , Subrahmanyam  (1971) . For 
more  complicated  problems,  it  is  necessary  to  resort  to  numerical 
methods;  the  use  of  the  finite  difference  solution  for  the  diffusion 
form  of  the  equations  of  consolidation  was  pioneered  by  Gibson  and 
Lumb  (1953) . 

The  diffusion  approach  is  only  an  approximation  to  the  actual 
process  of  consolidation,  because  it  neglects  the  changes  in  bulk  stress 
which  occur  during  consolidation.  % more  correct  formulation  for  a soil 
with  an  elastic  skeleton  was  developed  by  Biot  (1941a,  b,  1956a,  b) . 

Analytical  solutions  of  Biot's  equations  have  been  obtained  by  Biot 
(1955) (1956) , Josselin  De  Jong  (1957) (1963) , McNamee  and  Gibson  (1960a, 
b) , and  Gibson,  Schiffman  and  Pu  (1970)  for  particular  boundary^  condi- 
tions. However  such  solutions  are  limited  to  problems  with  extremely 
simple  geometry  and  boundary  conditions.  For  more  complicated  problems 
it  is  necessary  to  find  solutions  by  numerical  techniques.  Finite 
element  methods  have  been  developed  for  this  purpose  by  Sandhu  and  Wilson 
(1969),  Yokoo,  Yamagata  and  Nagaota  (1971),  Christian  and  Boehmer  (1970) 
and  Hwang  et  al.  (1971). 
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The  previous  theories  and  analyses  assume  that  the  soil  skeleton 
behaves  as  a linear  elastic  porous  medium.  Under  many  loading  conditions, 
however,  soils  do  not  behave  elastically.  Zienkiewicz  and  Naylor  (1972) , 
Small,  Booker  and  Davis  (1976)  have  developed  formulations  and  finite 
element  solution  techniques  for  the  consolidation  of  soils  with  elasto- 
plastic  stress-strain  behavior  of  the  soil  skeleton. 

For  analyses  of  consolidation  in  earth  dams,  it  is  important  to  be 
able  to  predict  the  pore  pressures  and  the  rate  at  which  the  pore  pressures 
dissipate  during  construction,  for  the  following  reasons: 

(1)  Considerable  movement,  both  vertical  and  horizontal,  may  occur 
due  to  pore  pressure  dissipation  and  may  cause  the  loss  of 
freeboard. 

(2)  Pore  pressure  generation  and  dissipation  leads  to  stress 
changes  during  consolidation  which  may  influence  the  likeli- 
hood of  hydraulic  fracturing  after  the  reservoir  is  filled. 

(3)  It  is  useful  to  be  able  to  predict  movements  and  pore  pres- 
sures in  dams  in  order  to  understand  the  significance  of 
measurements  made  during  and  after  construction. 

(4)  In  order  to  perform  effective  stress  analyses  of  stability, 
it  is  essential  that  the  pore  pressures  within  the  dam  be 
known. 

Previously  developed  procedures  for  analysis  of  consolidation  are 
only  applicable  to  saturated  clays.  Because  dams  generally  consist  of 
unsaturated  materials,  it  is  therefore  important  to  develop  a means  of 
describing  and  analyzing  the  consolidation  of  unsaturated  clay  soils. 

This  study  shows  that  the  air-water  mixture  in  the  pores  of  a clay 
compacted  wet  of  optimum  water  content  can  be  treated  as  a homogenized 
compressible  fluid.  Biot's  theory  can  be  used  to  analyze  this  condition 
if  the  following  factors  are  taken  into  account: 

(1)  the  effective  stress-strain  behavior  of  the  clay, 

(2)  the  compressibility  of  the  pore  fluid,  and 

(3)  the  permeability  of  the  unsaturated  clay. 
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These  factors  are  discussed  in  detail  in  this  study  and  are  used  to 
formulate  a theory  of  consolidation  applicable  to  the  behavior  of  unsatu- 
rated clays. 

The  final  part  of  this  study  is  concerned  with  the  development  of 
a numerical  technique  based  on  the  finite  element  method,  and  the  use  of 
this  technique  to  analyze  the  behavior  of  New  Melones  Dam. 
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CHAPTER  2 

EFFECTIVE  STRESS-STRAIN  RELATIONSHIP  OF  CLAY 

INTRODUCTION 

It  is  well  known  that  the  effective  stress  controls  the  failure  and 
deformation  behavior  of  soil,  and  so  it  is  important  to  investigate  the 
relation  between  effective  stress,  strains,  and  failure. 

The  first  attempt  to  describe  the  effective  stress-strain  relation- 
ship assumed  linear  elastic  behavior.  Although  this  approach  is  some- 
times applicable,  it  cannot  duplicate  the  irrecovable  deformation  observed 
in  real  soils,  nor  can  it  predict  the  possibility  of  failure. 

More  realistic  behavior  can  be  simulated  by  using  a non-linear 
elastic  model  with  a varying  Young's  modulus  and  Poisson's  ratio.  This 
type  of  model  predicts  real  soil  behavior  more  accurately  for  some 
conditions.  However  such  models  assume  that  the  principal  directions  of 
the  incremental  stress  and  incremental  strains  coincide,  which  is  general- 
ly incorrect  near  failure. 

One  way  to  overcome  these  limitations  is  to  use  the  theory  of 
plasticity.  A first  attempt  would  be  to  assume  that  the  soil  is  an 
elastic-perfectly  plastic  material.  Such  a model,  however,  does  not  pre- 
dict the  early  irrecovable  deformation  which  is  observed  in  practice. 
Drucker,  Gibson  and  Henkel  (1957)  recognized  that  in  order  to  obtain  a 
more  realistic  soil  behavior  it  was  necessary  to  propose  an  elastic- 
plastic  work  hardening  model. 

The  concept  of  soil  as  an  elastic-plastic  work  hardening  material 
has  been  studied  exhaustively  by  Roscoe  and  his  co-workers  in  Cambridge. 
The  Cam  Clay  Model  was  developed  by  Roscoe  and  Schofield  (1963)  and  Roscoe, 
Schofield  and  Thurairajah  (1963) . It  is  assumed  that  a conical  failure 
surface  is  capped  by  a yield  surface,  as  shown  in  Fig.  2.1,  and  that  the 
cap  yield  surface  expands  according  to  a work  hardening  law.  This 
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assumption  models  volume  changes  due  to  changes  in  both  shear  and  normal 
stress.  The  plastic  incremental  strain  is  assumed  to  be  governed  by  an 
associated  flow  rule.  The  shape  of  the  cap  yield  surface  was  later  modi- 
fied and  revised  by  Roscoe  and  Burland  (1968) . Similar  models  have  been 
investigated  by  Jenike  and  Shield  (1959),  Christian  (1966),  DiMaggio  and 
Sandler  (1971)  , and  Prevost  and  Hoeg  (1975) . 

The  main  differences  between  these  models  is  the  shape  of  the  caps. 
The  different  caps  which  have  been  used  are  summarized  below: 


(1) 

Spherical  cap:  Drucker,  Gibson  and  Henkel 

and  Shield  (1959) 

(1958) 

, Jenike 

(2) 

Bullet-shaped  cap:  Roscoe  and  Schofield  (1963) , 

Schofield  and  Thurairajah  (1963) 

Roscoe, 

(3) 

Elliptical  cap:  Christian  (1966) , Burland 

(1971) 

(1968) 

, Hagmamn 

(4) 

Exponential  cap:  DiMaggio  et  al.  (1971) 

(5) 

Conical  cap:  Prevost  and  Hoeg  (1975) . 

CAM  CLAY  MODELS 

This  study  is  concerned  with  the  stress-strain  behavior  of  com- 
pacted clays  in  the  core  of  an  earth  or  rockfill  dam.  The  Cam  Clay 
Models  were  developed  for  remolded  clays  and  have  been  used  successfully 
to  predict  their  stress-strain  behavior.  In  general,  the  properties  of 
compacted  clays  are  similar  to  those  of  remolded  clays,  and  it  thus  seems 
reasonable  to  adopt  the  concepts  of  the  Cam  Clay  Models  to  describe  the 
behavior  of  compacted  clays.  These  concepts  are  outlined  briefly  in  the 
following  section. 

Basic  Assumptions  on  Soil  Properties 

The  Cam  Clay  Model  assumes  that  the  failure  surface  can  be  expressed 
by  a linear  relation  between  a ' and  q: 
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(2.1) 


a ’ = Mq 
m 


where 


O'  = <ai  + + a,)/3 

m l z j 


i [Br« 


2)!+(02-o/  + (a3 


- V] 


The  failure  criterion  represented  by  Eq.  2.1  is  shown  in  Fig.  2.2. 

In  describing  the  volume  change  characteristics,  it  is  assumed  that 

the  isotropic  compression  curve  (the  e-  In  O ' curve)  can  be  represented 

m 

by  two  lines  with  different  slopes,  X and  K,  as  shown  in  Fig.  2.3. 

Elastic  Behavior 

An  elastic-plastic  model  behaves  elastically  during  unloading. 
Therefore  the  rebound  curve  of  an  isotropic  compression  curve  can  be  used 
to  determine  the  elastic  volume  change  behavior. 

The  relationship  between  an  incremental  change  in  void  ratio  (de) 
and  the  corresponding  increment  of  volumetric  strain  (dev)  is: 

dev = rrr;  (2-4) 


where  eQ  is  the  initial  void  ratio. 
Bulk  modulus  is  defined  as: 


With  the  aid  of  Eq.  2.4,  the  bulk  modulus  can  be  expressed  as: 


"de~  (1  + e0) 


As  shown  in  Fig.  2.3,  < is  the  slope  of  the  rebound  curve  in  the 

e - In  o ' plot.  That  is: 
m 
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1FIG.  2.3  SIMPLIFIED  e - In  o ' CURVE 

m 
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K 


(2.7) 


Substituting  Eq.  2.7  into  Eq.  2.6,  the  bulk  modulus  can  be  calcu- 
lated for  different  values  of  mean  effective  stress  a ' as  follows: 

m 


B = — (l  + e0)  (2.8) 

In  order  to  determine  the  elastic  behavior,  it  is  necessary  to 
determine  two  elastic  parameters.  For  example,  the  bulk  modulus  and  the 
shear  modulus  may  be  used.  In  order  to  obtain  an  analytically  simple 
formulation,  elastic  shear  strains  are  generally  neglected,  which  is 
equivalent  to  assuming  that  the  shear  modulus  is  infinite.  This 
assumption  is  not  strictly  necessary  as  it  is  possible  to  formulate  the 
theory  for  any  value  of  shear  modulus.  When  this  is  done,  it  is  observed 
that  the  stress-strain  behavior  is  not  particularly  sensitive  to  the  value 
of  shear  modulus.  In  many  situations,  it  is  permissible  to  select  a 
reasonable  value  of  Poisson's  ratio  in  the  range  0.2  £ V'  £ 0.45,  and  to 
use  the  corresponding  values  of  shear  modulus  in  the  analyses. 

Plastic  Behavior 

The  original  Cam  Clay  Model  developed  in  1963  has  a bullet-shaped 
yield  surface  which  can  be  expressed  as: 


<V  “*(3)  - Po(CvP)-  0 


where 


(2.10) 


pQ(evP\  = the  hydrostatic  pressure  necessary  to  produce 

P 

a plastic  volumetric  strain  of  magnitude  e . 


The  yield  surface  expressed  in  Eq.  2.9  expands  according  to  the 
work  hardening  parameter  P0(evP)'  whic^  can  be  obtained  by  deducting  the 
elastic  volumetric  strain  from  the  total  volumetric  strain.  Thus  the 
plastic  volumetric  strain  can  be  expressed  as 
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(2.11) 


a ' 

(A-  k)  In  , (a  • > a •) 

a ' m mp 

mp 


p 

where  0 ' is  the  value  of  mean  effective  stress  for  which  £ is  zero, 

mp  v 

as  shown  in  Fig.  2.3. 

In  the  Cambridge  theory  it  is  assumed  that  the  plastic  rtrain 
increment  is  given  by  an  associated  flow  rule.  Thus  tne  vector  of  plastic 
strain  increment  can  be  determined  by  the  shape  of  the  yield  cap.  For 
some  soils,  Roscoe  and  Burland  (1968)  suggested  that  an  elliptical  yield 
surface  may  be  more  appropriate  than  the  bullet  shape  yield  surface  in 
Eq.  2.9.  The  elliptical  yield  surface  can  be  expressed  as 

”2  +~2  ■ O'  - P (e  p)  = 0 (2.12) 

M2  m 0 v 

The  modified  Cam  Clay  theory  can  be  obtained  by  replacing  Eq.  2.9 
by  this  equation.  Both  yield  surfaces  are  shown  schematically  in  Fig. 

2.4. 


STRESS-STRAIN  BEHAVIOR  OF  COMPACTED  CLAY  COMPARED  WITH  THAT  PREDICTED  BY 
THE  CAM  CLAY  MODELS 

This  study  is  concerned  with  the  pore  pressure  dissipation  and  the 
stress-strain  behavior  of  compacted  clay  in  the  core  of  earth  and  rock- 
fill  dams.  Therefore  it  is  interesting  to  see  how  well  the  Cam  Clay 
Models  can  predict  the  observed  stress-strain  behavior  of  compacted  clays. 
In  order  to  do  this,  the  results  of  predictions  made  using  both  Cam  Clay 
Models  are  compared  with  the  observed  behavior  of  compacted  Canyon  Dam 
clay  which  was  tested  by  Casagrande  (1962) . 

Canyon  Dam  clay  is  a silty  clay  with  PI  = 19  and  LL  = 34.  The 
maximum  dry  density  as  determined  by  the  Harvard  compaction  test  is 
112.8  lb/ft3  and  the  optimum  water  content  is  16.7%.  Test  samples  of 
this  clay  were  compacted  2%  wet  of  optimum  water  content  with  an  average 
water  content  17.4%,  and  an  average  dry  density  111.3  lb/ft3. 

The  tests  available  for  parameters  evaluation  were: 
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(1)  Drained  triaxial  tests,  and 

(2)  One-dimensional  compression  tests. 

The  results  of  these  tests  are  shown  in  Fig.  2.5  and  Fig.  2.6. 

The  value  of  the  parameter  M can  be  determined  from  the  drained 
triaxial  tests.  The  failure  line  from  the  drained  triaxial  tests  is 
plotted  in  Fig.  2.7. 

It  is  noted  that  the  most  appropriate  failure  line  for  this  clay 
does  not  pass  through  the  origin  as  postulated  in  the  Cam  Clay  Models. 

It  is  most  desirable  to  determine  X and  K from  an  isotropic  com- 
pression curve.  Unfortunately,  this  type  of  test  was  not  available  for 
Canyon  Dam  clay.  However,  the  results  of  a one-dimensional  compression 
curve  were  available.  Therefore,  the  assumption  that  KQ  remained  con- 
stant through  the  test  was  made  so  that  the  slopes  X and  K measured  from 
the  one-dimensional  compression  test  could  be  assumed  to  be  the  same  as 
those  measured  from  isotropic  consolidation  tests.  Values  of  Poisson’s 
ratio  could  be  determined  if  the  results  of  drained  triaxial  tests  on 
overconsolidated  specimens  were  available.  Such  tests  had  not  been  con- 
ducted for  Canyon  Dam  clay,  but,  as  observed  before,  the  stress - 
strain  behavior  is  not  sensitive  to  the  value  of  Poisson's  ratio.  Thus 
a value  of  0.3  was  chosen  for  V'  and  was  used  in  the  calculations. 

By  using  the  parameters  summarized  in  Table  2.1,  the  stress-strain 
curves  for  the  drained  triaxial  tests  cam  be  predicted.  Comparisons 
between  the  predicted  and  the  experimental  stress-strain  curves  are  shown 
in  Fig.  2.8.  It  may  be  seen  that  the  stress-strain  behavior  predicted  by 
the  elliptical  cap  (Roscoe  and  Burland,  1968)  is  somewhat  better  than 
that  predicted  by  the  bullet  shaped  cap  (Roscoe  and  Schofield,  1963) . 
However,  the  comparison  shows  the  stress-strain  behavior  predicted  by 
these  theories  is  not  altogether  satisfactory. 

The  predicted  stress-strain  behavior  is  deficient  in  two  ways: 
the  strength  is  not  predicted  accurately,  and  the  calculated  values  of 
shear  strain  C are  too  large. 
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Axiol  Strain,  C(  (%) 


FIG.  2.5  STRESS-STRAIN  CURVES  FROM  DRAINED 
TRIAXIAL  TESTS  ON  CANYON  DAM 
COMPACTED  CLAY 

(After  Casagrande,  1962) 
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FIG.  2.6  ONE  DIMENSIONAL  CONSOLIDATION  CURVE  FOR  CANYON  DAM  COMPACTED  CLAY 

(After  Casagrande,  1962) 


FIG.  2.7  FAILURE  LINES  FOR  CANYON  DAM  COMPACTED  CLAY 
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Table  2.1  Summary  of  the  Parameters  for 
Cam  Clay  Models  (1963  & 1968) 


Parameters 

Cam  Clay  (1963) 

Cam  Clay  (1968) 

X 

0.05 

0.05 

K 

0.004 

0.004 

M 

1.25 

1.25 

V' 

0.3 

0.3 

17 


FIG.  2.8  COMPARISONS  OF  MEASURED  AND  CALCULATED 
STRESS-STRAIN  CURVES 
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The  possible  reasons  for  this  are: 

(1)  The  strength  of  Canyon  Dam  clay  is  not  well  represented  by  a 
failure  line  passing  through  origin,  as  shown  in  Fig.  2.7. 

(2)  The  use  of  two  straight  lines  to  simulate  the  measured 

e - In  0 ' curve  may  not  be  sufficiently  accurate, 
m 

The  following  section  describes  a modification  of  the  Cam  Clay  Model 
intended  to  overcome  these  defects. 

A MODIFICATION  OF  THE  CAM  CLAY  MODEL 

In  order  to  be  able  to  predict  more  realistic  stress-strain  beha- 
vior for  compacted  clays,  the  following  modifications  were  made  to  the 
Cam  Clay  Model. 

(1)  It  was  assumed  that  the  failure  line  had  a cohesion  intercept 

so  that  the  strength  of  the  soil  could  be  represented  more 

accurately.  The  failure  line  in  the  0 ' - q plane  may  be 

m 

expressed  as 

t 

q = qQ  + 0^'  M (2.13) 

where  qQ  is  the  intercept  on  the  q axis  of  the  q - 0^'  diagram. 

(2)  The  measured  e - In  O ' curve  is  used  instead  of  a straight 

in 

line  approximation.  The  e - In  o^'  relationship  was  repre- 
sented by  a table  of  pairs  of  values  of  e and  o ' , and  inter- 
mediate values  were  determined  by  interpolation. 

Elastic  Behavior 

With  a rebound  curve  represented  by  pairs  of  values  of  e and  o • , 

m 

the  tangent  value  of  bulk  modulus  can  be  expressed  as 

B - (1  + eQ)/b  (2.14) 

where  b is  the  tangent  value  of  the  slope  of  the  rebound  e-0  ' curve. 

m 

For  Canyon  Dam  compacted  clay  the  value  of  B is  very  nearly  constant 
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over  a range  of  values  of  o'  as  shown  in  Fig.  2.9.  The  value  of 
Poisson's  ratio  was  assumed  to  be  0.3,  as  discussed  previously. 


Plastic  Behavior 

With  a failure  line  which  does  not  pass  through  the  origin,  the 
elliptical  yield  cap  can  be  expressed  as 

(“)(^  r)--r-b(SP) 

where  p^  = q^/M  (2.16) 

p 0(evP)  can  be  obtained  from  a o^'  - e curve  by  deducting  the 
elastic  volumetric  strain  from  the  total  volumetric  strain.  The  plastic 
strain  increments  are  calculated  using  an  associated  flow  rule. 


Parameters  Required 


The  parameters  required  to  implement  the  theory  are  B,  V' , M,  c and 


the  pairs  of  values  of  e and  0 ' which  represent  the  isotropic  compres- 

m 

sion  curve.  Values  of  these  parameters  for  Canyon  Dam  clay  are  evalu- 
ated in  the  following  section. 


STRESS-STRAIN  BEHAVIOR  OF  CANYON  DAM  CLAY  PREDICTED  USING  THE 
EXTENDED  CAM  CLAY  MODEL 

The  values  of  the  parameters  M and  qQ  can  be  determined  from 
Fig.  2.7.  The  isotropic  compression  curve  was  not  available,  and  there- 
fore, as  before,  the  0 ' - e curve  was  calculated  using  the  one-dimensional 

m 

compression  curve  shown  in  Fig.  2.6.  In  calculating  values  of  e and  0^' 
for  the  curve  shown  in  Fig.  2.9,  it  was  assumed  that  KQ  was  equal  to  0.5. 

The  evaluation  of  the  slope  of  the  rebound  curve  (b)  is  shown  in 
Fig.  2.9. 

Using  the  parameters  summarized  in  Table  2.2,  the  Extended  Cairn  Clay 
Model  was  used  to  predict  the  stress-strain  curves  for  Canyon  Dam  clay. 
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0.5 


0 5 10  15 

Isotropic  Pressure,  (kg/cm2) 


FIG.  2.9  ISOTROPIC  COMPRESSION  CURVE 
FOR  CANYON  DAM  CLAY 
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Table  2.2  Summary  of  the  Parameters  for 
the  Extended  Cam  Clay  Model 


Parameters 

Values 

B 

2250  kg/cm3 

V' 

0.3 

M 

1.18 

c 

0.9  kg/cm3 

p0(e) 

in  Fig.  2.9 
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As  shown  in  Fig.  2.10,  the  Extended  Model  provides  significant  improve- 
ment in  the  prediction  of  both  strength  and  the  shape  of  the  stress-strain 
curves. 
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MODIFICATIONS  TO  THE  SHAPE  OF  THE  CAP 

The  plastic  strain  increments  are  determined  by  the  shape  of  the 
cap.  For  some  soils,  the  stress-strain  behavior  cannot  be  described 
sufficiently  well  by  the  models  described  previously.  Therefore,  some 
further  modifications  were  considered  to  extend  the  range  of  applicability 
of  the  elastic-plastic  work  hardening  model. 

Curved  Failure  Envelope  for  Granular  Soils 

For  some  granular  materials,  it  is  found  that  the  failure  envelope 
is  curved  due  to  particle  crushing  at  high  confining  pressures.  The 
corresponding  reductions  in  the  values  of  <t>'  with  increasing  pressure  may 
be  represented  by  the  equation  proposed  by  Wong  and  Duncan  (1974) : 

°3 

4>*  = 4>0'  - A<f>'  iog1Q  — (2.17) 

pa 

where  p is  the  atmospheric  pressure,  <J>  • is  the  value  of  <j>'  for  O' 
equal  to  p , and  A<f>’  is  the  reduction  in  <p'  for  a ten-fold  increase  in 

3 

03' . Eq.  2.17  can  be  used  to  evaluate  the  angle  of  internal  friction 
appropriate  to  any  confining  pressure  within  the  range  of  pressures 
encompassed  by  the  test  results.  The  slope  of  failure  line  in  o^-  q 
plane,  M,  can  be  related  to  by  the  equation: 


„ 6sinJ' 

3 - sin  <P' 


(2.18) 


Therefore  the  cap  is  a function  of  M,  and  the  value  of  M varies 
with  confining  pressure  o^  . 
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Axiol  Strain,  £,  (%) 


FIG.  2.10  COMPARISON  OF  MEASURED  STRESS-STRAIN 
CURVES  WITH  THOSE  CALCULATED  USING 
THE  EXTENDED  CAM  CLAY  MODEL 


Composite  Cap 


For  some  soils  the  observed  shear  strain  is  larger  than  that  pre- 
dicted by  the  Cam  Clay  Models  discussed  previously.  This  limitation  can 
be  overcome  by  a proposed  revision  of  the  cap  suggested  by  Roscoe  and 
Burland  (1968). 


According  to  the  concept  proposed  by  Roscoe  and  Bur land,  it  is 
assumed  that  the  yield  surface  consists  of  two  parts,  AB  and  BC,  as 
shown  in  Fig.  2.11.  The  portion  AB  is  elliptical,  and  BC  is  a straight 
line.  The  elliptical  part  can  be  expressed  by  Eq.  2.12.  Another  harden- 
ing parameter  E13  is  chosen  to  describe  the  straight  line  yield  surface. 
e*3  is  the  octahedral  plastic  shear  strain,  which  is  given  by  the  follow- 
ing equation: 


■ i y 2 [(eip  - 


e2P) 2 + (e2P  - e3P)  2 + (e3P  - e^)  2 


(2.19) 


Experimental  results  suggest  that  the  yield  function  of  the  straight 
line  part  of  the  cap  surface  can  be  expressed  as: 


a 


(Mp)  1 


- £p  = o 


(2.20) 


where  a is  a constant. 

It  is  assumed  that  an  associated  flow  rule  can  be  used  for  both 
parts  of  the  cap  yield  surface  for  the  determination  of  the  plastic  strain 
increment.  On  the  portion  BC,  only  plastic  shear  strain  occurs,  without 
any  plastic  volume  strain. 

At  the  corner  B,  where  the  two  parts  of  the  yield  surface  intersect, 
it  is  assumed  that  the  plastic  strain  increment  can  be  obtained  by  super- 
imposing plastic  strains  derived  from  the  two  yield  surfaces  AB  and  BC 
(Koiter,  1953).  The  plastic  strain  increments  can  be  expressed  as: 

d£  P ■ d£  P (2.21) 

V V 

AB 

d?  - de*  + d?3  (2.22) 

AB  BC 
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in  which  the  de  P and  de*3  are  determined  from  the  elliptical  yield 
v.  D AB 

AB  ~p 

surface  AB  , and  de  is  determined  from  the  yield  surface  BC. 

BC 

The  parameter  a governs  the  magnitude  of  de*3  . One 

BC 

possible  way  of  determining  the  parameter  a is  to  determine  values  of 

de*3^  by  subtracting  the  corresponding  values  of  de*3^  and  the  elastic 

strain  increments  from  the  total  strains.  Then  the  appropriate  a could 

be  calculated  using  Eq.  2.20.  This  procedure  is  not  simple,  however, 

because  the  values  of  de*3,..  are  difficult  to  calculate. 

AB 

Perhaps  the  most  straight-forward  way  of  determining  a for  practi- 
cal purposes  is  by  curve  fitting.  If  all  the  other  parameters  are  known, 
the  stress-strain  curves  for  a series  of  drained  triaxial  tests  can  be 
calculated  for  different  values  of  a,  and  a value  of  a which  best  fits 
the  tests  can  be  chosen.  A computer  program  for  determining  the  stress- 
strain  curves  for  triaxial  tests  is  provided  in  Appendix  C.  An  example 
is  shown  in  Appendix  D. 
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CHAPTER  3 


THE  PERMEABILITY  OF  PARTLY  SATURATED  SOILS 


INTRODUCTION 

The  permeability  of  the  compacted  clay  core  of  an  earth  dam 
depends  on  a number  of  factors,  including  the  type  of  soil,  the  compac- 
tion conditions,  the  void  ratio  or  density,  and  the  degree  of  saturation. 
The  coefficient  of  permeability  may  be  measured  in  laboratory  tests  on 
samples  contacted  under  the  same  conditions  as  those  in  the  field.  How- 
ever, the  permeability  will  change  as  the  degree  of  saturation  changes 
during  seepage,  and  it  will  also  change  as  the  clay  is  compressed  under 
the  weight  of  overlying  fill.  In  an  analysis  of  the  consolidation  of 
the  core  of  a dam,  it  is  important  to  take  into  account  the  variations 
of  permeability  with  void  ratio  and  degree  of  saturation. 

In  this  chapter,  an  empirical  equation  is  proposed  for  describing 
the  changes  in  permeability  of  partly  saturated  compacted  clays  which 
result  from  changes  in  degree  of  saturation  and  void  ratio. 


THEORETICAL  PERMEABILITY  EQUATIONS 

There  have  been  many  attempts  at  establishing  relationships  between 
soil  permeability  and  other  soil  properties.  Poiseiulle  developed  a 
formula  for  the  permeability  of  a porous  medium  with  irregular  pores.  He 
found  that 


(3.1) 


where  k is  the  coefficient  of  permeability,  p is  the  viscosity  of  fluid, 
is  the  unit  weight  of  water  (these  are  both  functions  of  temperature) , 
Vs  is  the  volume  of  solid,  Ag  is  the  wetted  area,  e is  the  void  ratio, 

S is  the  degree  of  saturation  and  is  a shape  coefficient  whose  value 
depends  on  the  degree  of  irregularity  of  the  pore  spaces. 
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Formulae  similar  to  Poiseiulle's  equation  have  been  also  developed 
by  Kozney  (1927),  Taylor  (1948),  and  Carman  (1956).  All  of  these  theo- 
retical equations  are  based  on  the  assumptions  that  the  particles  are 
approximately  equidimensional , that  they  are  of  a uniform  size  larger 
than  1 urn,  and  that  the  flow  is  laminar.  The  equations  have  been 
observed  to  hold  fairly  well  in  sandy  and  silty  soils.  However,  consider- 
able discrepancies  between  these  theoretical  expressions  and  experimental 
results  have  been  found  for  clays.  It  seems  well  established  that  the 
major  cause  of  these  discrepancies  is  that  the  pores  of  clay  soils  are  not 
of  a uniform  size.  Olsen  (1961)  has  discussed  these  discrepancies  in 
detail. 


PROPOSED  PERMEABILITY  EQUATION 

Because  the  theoretical  permeability  equations  are  not  applicable 
to  clays,  it  is  preferrable  to  adopt  an  experimental  or  empirical 
approach . 

Soil  permeability  is  influenced  by  many  factors  including  particle 
size  distribution,  particle  shape,  structure,  void  ratio  and  degree  of 
saturation.  The  effects  of  each  of  these  is  hard  to  isolate  from  the 
others  because  these  characteristics  are  closely  interrelated — e.g., 
structure  depends  on  particle  size,  void  ratio,  method  of  compaction, 
and  compaction  water  content. 

In  the  study  described  in  this  report  no  attempt  has  been  made  to 
formulate  the  effects  of  particle  size  distribution,  particle  shape  and 
structure  on  permeability,  because  these  factors  are  nearly  constant  for 
a particular  clay,  a particular  method  of  compaction  and  a given  compac- 
tion water  content,  and  their  effects  are  accounted  for  when  laboratory 
tests  are  performed  on  specimens  compacted  under  the  same  conditions  as 
the  soil  in  the  field. 

Two  major  factors  which  change  the  permeability  of  compacted  core 
material  after  compaction  are  changes  in  the  degree  of  saturation  and  the 
void  ratio.  It  seems  reasonable  to  postulate  that  their  effects  can  be 
accounted  for  by  using  an  empirical  equation  of  the  form: 
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k = k G H (3.2) 

s e s 

in  which  k is  the  permeability  of  the  partly  saturated  clay,  kg  is  the 
permeability  of  the  clay  when  it  has  been  completely  saturated  using 
back  pressure,  Ge  is  a factor  whose  value  depends  on  void  ratio,  and  Hg 
is  a factor  whose  value  depends  on  degree  of  saturation. 

Determination  of  k 
s 

The  value  of  kg  is  influenced  significantly  by  the  water  content 
at  compaction  because  the  water  content  at  compaction  influences  the 
structural  arrangement  of  soil  properties.  Therefore  kg  must  be  deter- 
mined from  test  specimens  compacted  using  procedures  which  produce  the 
same  structure  as  does  the  field  compaction  procedure.  Seed  and  Chan 
(1959)  found  that  a flocculated  structure  ( random  particle  arrangement) 
tends  to  be  formed  for  clays  compacted  dry  of  optimum,  while  a dispersed 
structure  (more  parallel  particle  arrangement)  tends  to  be  formed  for 
clays  compacted  wet  of  optimum  by  kneading  compaction.  As  shown  in 
Fig.  3.1,  the  permeability  of  a soil  with  a flocculated  structure  is 
far  higher  than  the  permeability  of  the  same  soil  with  a dispersed 
structure . 


Determination  of  G 
e 

The  variation  of  the  factor  G with  void  ratio  can  be  determined 

e 

experimentally  by  testing  saturated  specimens  which  have  been  consoli- 
dated under  different  pressures,  and  therefore  have  different  void  ratios. 


Data  from  a number  of  tests  on  a sand  are  shown  in  Fig.  3.2.  For 
this  sand,  the  variation  of  kg  with  e appears  to  be  fairly  accurately 
approximately  by  Poiseiulle's  equation  which  corresponds  to 


G 

e 


e3  / eQ3 
1 + e / 1 + eQ 


(3.3) 


Additional  data  for  a number  of  sands,  silts,  and  clays  are  shown 
in  Fig.  3.3  (from  Lambe  and  Whitman,  1969).  These  data  indicate  that 
the  log  of  kg  varies  linearly  with  e.  This  variation  corresponds  to 
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FIG.  3.1  PERMEABILITY  BEHAVIOR  OF  JAMAICA  SANDY  CLAY 

(After  Lanbe,  1958) 


3 


Void  Rotio  Function 


FIG.  3.2  VARIATION  OF  PERMEABILITY  WITH 
VOID  RATIO  FROM  PERMEABILITY 
TESTS  ON  SAND 

(After  Lambe  and  Whitman,  1969) 
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PERMEABILITY  TEST  DATA  (After  Lambe  and  Whitman,  1969) 


(3.4) 


Ge  = Exp(|)/  ExP  (T) 

where  8 is  the  slope  of  the  e-  log  kg  line. 

For  practical  purposes,  either  Eq.  3.3  or  Eq.  3.4  could  be  used  to 
define  Ge»  depending  on  which  best  represents  the  observed  behavior  of 
the  particular  compacted  clay  under  consideration.  If  experimental  data 
are  not  available  for  the  soil  under  consideration,  it  seems  reasonable 
to  assume  that  Ge  is  defined  by  Eq.  3.3. 

Determination  of  H 
s 

Experimental  data  (Singh,  1965)  suggest  that  the  relationship 
between  the  permeability  of  clays  and  their  degree  of  saturation  can 
often  be  approximated  using  this  following  relationship  between  Hg  and 
S. 


where  Sf  is  the  threshold  degree  of  saturation  at  which  the  water  in  the 

pores  begins  to  flow  freely.  For  a degree  of  saturation  less  than  S^, 

the  water  forms  a membrane  around  the  particles  and  does  not  flow. 

When  m = 3 and  S.  = 0,  Eq.  3.5  conforms  to  the  variation  of  H with  S 
t s 

implicit  in  Poiseiulle's  equation.  The  value  of  varies  depending  on 
the  structure  of  the  clay.  Therefore  compaction  water  content  and 
contraction  method  have  an  influence  on  the  value  of  Sf. 

The  curves  plotted  in  Fig.  3.4  indicate  that  values  of  Sf  in  the 
range  from  Sf  * 0.0  to  = 0.6,  together  with  m = 3,  approximate  closely 
the  results  of  tests  conducted  by  Mitchell  et  al.  (1965) . 

If  experimental  data  are  not  available  to  define  the  most  appro- 
priate value  of  H , it  seems  reasonable  to  assume  that  Poiseiulle's 
s 

equation  is  applicable,  i.e. , that  Sf  * 0 and  m = 3.  If  this  is  done, 

Hg  will  be  related  to  S by  the  following  equation: 

Hg  - S3  (3.6) 
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FIG 


. 3.4  VARIATIONS  OF  PERMEABILITY  WITH  DEGREE 
OF  SATURATION  FOR  COMPACTED  CLAY 

(Data  from  Mitchell,  et  al.,  1965 


SUMMARY 


The  permeability  of  compacted  clay  depends  on  the  type  of  soil, 
the  compaction  conditions,  the  degree  of  saturation,  and  the  void  ratio. 
The  effects  of  changes  in  degree  of  saturation  and  void  ratio  can  be 
approximated  using  Eq.  3.2.  The  value  of  saturated  permeability,  kg,  can 
be  measured  using  test  specimens  compacted  under  the  same  condition  as 
those  in  the  field. 

The  variation  of  Gg  with  void  ratio  can  be  determined  experimental- 
ly provided  the  necessary  test  results  are  available.  If 
available,  it  appears  to  be  reasonable  to  assume  that  Gg  = 

as  suggested  by  Poseiulle's  equation.  This  is  what  was  done  for  the 
analysis  of  New  Melones  Dam. 

The  variation  of  Hg  with  S can  also  be  determined  ejqjerimentally 
if  test  results  are  available.  If  no  data  are  available  to  define  this 
variation,  it  can  be  assumed  again  that  Poiseiulle's  equation  is  appli- 
cable. This  assumption  corresponds  to  H = S3,  which  has  been  found  to 

s 

agree  reasonably  well  with  experimental  data  obtained  by  Mitchell  et  al., 
(1965) . Because  no  data  were  available  to  define  the  variation  of  H 

s 

with  S for  the  New  Melones  core  material,  it  was  assumed  that  Hg=S3. 


data 

03 


are  not 


0 


1 + e 
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CHAPTER  4 

EFFECTIVE  STRESS  EQUATION  FOR  PARTLY  SATURATED  SOILS 
INTRODUCTION 

It  is  well  established  that  effective  stresses  govern  the  defor- 
mation and  strength  of  soils.  For  a saturated  soil,  the  expression  for 
effective  stress  may  be  written  as: 

O'  = O - u (4.1) 

where  a is  total  stress,  u is  the  pore  fluid  pressure  and  O'  denotes 
effective  stress. 

For  an  unsaturated  soil,  the  pore  fluid  consists  of  air  and  water. 
The  effective  stress  equation  involves  water  pressure  u^,  air  pressure 
u , and  surface  tension  between  air  and  water.  Considering  equilibrium 
of  force  in  a microscopic  section  of  soil , the  effective  stress  equation 
can  be  expressed  as: 

The  total  forces  = the  interparticle  forces 

+ the  sum  of  the  water  pressure  forces 

+ the  sum  of  the  air  pressure  forces 

+ the  surface  tension  forces  (4.2) 

The  effective  stress  equation  is  thus  quite  simply  conceptually. 
However,  it  is  difficult  to  express  quantitatively. 

A possible  form  of  the  effective  stress  equation  for  partly  satu- 
rated soils  had  been  suggested  by  Bishop  (1959) : 

0'  = 0 - u +X(u  - u ) (4.3) 

a a w 

where  x is  an  empirical  parameter  representing  the  portion  of  the  surface 

tension,  (u  - u ) , that  contributes  to  the  effective  stress, 
a w 

X is  a function  of  degree  of  saturation.  The  form  of  this  relation 
has  been  derived  by  Aitchison  (1960) , who  assumed  that  the  soil  consists 
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of  uniformly  packed  spheres.  This  relation  is  plotted  in  Fig.  4.1.  How- 
ever, a number  of  discrepancies  have  been  found  between  the  theoretical 
curve  and  experimental  data  for  clays  (Blight,  1965) . Also  it  has  been 
found  that  the  values  of  X measured  from  strength  behavior  may  be  signifi- 
cantly different  from  those  measured  from  volume  change  behavior  (Coleman, 
1962;  Matyas  and  Radhakrishna,  1968) . 

Another  possible  relationship,  suggested  by  Sparks  (1963),  is: 

0'  = 0 - ft  u - £ u + ipT  (4.4) 

awe 

where  Tcis  the  surface  tension  and  ft-  £ and  are  parameters  represent- 
ing the  portions  of  u , u , and  T which  contribute  to  the  effective 

cl  W C 

stress. 

Sparks  employed  a model  composed  of  uniform  spheres  arranged  in 
various  packings,  and  having  various  contact  angles  between  the  water 
and  the  soil  particles.  On  the  basis  of  these  studies.  Sparks  concluded 
that  X may  assume  values  greater  them  one  or  less  than  zero,  as  shown  in 
Fig.  4.1 

Sparks  also  developed  an  extensive  set  of  charts  to  be  used  to 
determine  the  values  ft,  £ and  for  different  packing  conditions  and 
contact  angles.  All  of  the  values  proposed  by  Sparks  were  calculated  on 
the  assumption  that  the  soil  consists  of  uniformly  packed  spheres. 
Therefore  the  values  of  these  parameters  are  more  accurate  for  sandy  and 
silty  soils  than  for  clayey  soils. 

A frequent  special  case  occurs  when  the  clay  is  compacted  at  high 
water  contents,  and  occluded  air  bubbles  form  in  the  pores.  Under  this 
condition,  Sparks  assumed 

O'  = a - (1  - A)  u (4.5) 

w 

where  A is  the  effective  soil  particle  contact  area  per  unit  area.  At 
high  water  contents,  the  area  of  contact  between  the  solid  particles  is 
small  compared  to  the  area  of  contact  between  the  solid  particles  and  the 
water.  Under  this  condition,  the  value  of  A approaches  zero.  With  A=0, 
Eq.  4.5  becomes 
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FIG.  4.1  VARIATIONS  OF  X WITH  DEGREE 
OF  SATURATION 
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If 

r 

a'  = a - u (4.6) 

w 

Thus,  according  to  Sparks'  studies,  the  same  effective  stress 
equation  which  applies  to  saturated  soils  can  also  be  used  for  partly 
saturated  conditions  where  the  degree  of  saturation  is  high  and  the  air 
in  the  voids  is  in  the  form  of  occluded  bubbles. 

CONCEPT  OF  HOMOGENIZED  PORE  FLUID 

During  the  investigation  described  in  this  study,  it  was  concluded 
that  the  effective  stress  equation  previously  developed  for  partly  satu- 
rated soils  were  either  inaccurate  or  overly  complex.  In  an  attempt  to 
find  an  alternative,  the  concept  of  an  "homogenized  pore  fluid"  was 
developed. 

The  concept  is  that,  although  the  pore  fluid  of  a partly  saturated 
material  consists  of  a gas  and  a liquid,  it  can  be  treated  phenomenologi- 
cally as  a single  homogenized  pore  fluid.  The  pore  spaces  are  viewed  as 
being  completely  filled  with  this  homogenized  fluid,  and  the  properties 
of  the  fluid  change  as  the  portions  of  air  and  water  vary.  The  effective 
stress  law  corresponding  to  this  concept  of  an  homogenized  pore  fluid  is 
expressed  as: 

o'  = a - u (4.7) 

m 

where  u^  is  the  effective  average  pressure  in  the  homogenized  pore  fluid. 

For  an  occluded  bubble  system,  u^  is  simply  u^,  as  shown  by  Sparks.  For 

clays  with  lower  degrees  of  saturation,  u is  a function  of  the  air 

m 

pressure,  the  water  pressure  and  the  surface  tension. 

THE  HOMOGENIZED  PORE  FLUID  PRESSURE 

Because  the  three-phase  unsaturated  soil  is  treated  as  a saturated 

material  with  an  homogenized  pore  fluid,  the  total  stress  is  the  sum  of 

effective  stress.  O' , and  pressure  in  the  homogenized  pore  fluid,  u , as 

m 

expressed  by  Eq.  4.7. 

0=0'  + u (4.8) 

m 

i 

i 
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The  effective  stress  equation  for  partly  saturated  soil  suggested 
by  Bishop  (1959)  corresponds  to  the  following  homogenized  pore  fluid 
pressure : 

Um  = ua  ' X (ua  ~ Uw>  (4-9) 

The  effective  stress  equation  suggested  by  Sparks  (1963)  corres- 
ponds to  the  following  homogenized  pore  fluid  pressure: 

u = ft  u + £ u + i^T  (4.10) 

It  has  been  observed  previously  that  neither  of  there  effective 
stress  equations  is  entirely  satisfactory  for  clays  with  low  degrees  of 
saturation.  Fortunately,  the  clay  used  as  the  core  material  in  an  earth 
dam  is  usually  compacted  at  a water  content  corresponding  to  a degree  of 
saturation  in  the  neighborhood  of  85%.  In  this  case,  the  air  bubbles 
are  likely  to  be  occluded.  Therefore,  as  observed  by  Sparks,  it  can  be 
assumed  with  a high  degree  of  accuracy  that  u^  is  equal  to  u . This 
was  done  in  the  analysis  of  the  consolidation  of  the  core  of  New  Melones 
Dam. 


ANALYSIS  OF  CONSOLIDATION  OF  PARTLY  SATURATED  SOILS 

As  a partly  saturated  soil  is  loaded,  the  skeleton  is  subjected  to 
changes  in  the  effective  stress  O'  and  undergoes  deformations.  At  the 
same  time,  the  homogenized  pore  fluid  is  subjected  to  changes  in  pore 
pressure,  and  it  compresses  or  expands.  Under  the  influence  of  the 
hydraulic  gradient  induced,  the  pore  fluid  will  flow. 

The  homogenized  pore  fluid  model  can  be  used  for  analyses  of  con- 
solidation provided  the  following  are  known: 

(1)  The  relationship  between  effective  stress  and  strain  for  the 
soil  skeleton. 

(2)  The  relationship  between  permeability  and  the  degree  of  satu- 
ration of  the  homogenized  pore  fluid. 

(3)  The  relationship  between  the  compressibility  of  the  homogenized 
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pore  fluid  and  the  proportions  of  air  and  water  it  con- 
tains. 

The  effective  stress-strain  relationship  was  discussed  in  Chapter 
2,  and  the  permeability  relationship  was  discussed  in  Chapter  3.  The 
compressibility  relationship  will  be  discussed  in  the  following  sections 
of  this  chapter. 


THE  COMPRESSIBILITY  OF  THE  HOMOGENIZED  PORE  FLUID  OF  PARTLY 
SATURATED  SOIL 


The  compressibility  of  the  homogenized  pore  fluid,  ^ , is  defined 
as  the  volumetric  strain  induced  by  a unit  change  in  pore  fluid  pressure. 
It  ran  be  defined  by  the  following  equation: 


1 

Q 


de 

v 


(4.11) 


where  e is  the  volumetric  strain  and  u is  the  pore  fluid  pressure, 
v m 

As  observed  previously,  the  homogenized  pore  fluid  pressure  can  be 

expressed  as  a function  of  air  pressure,  water  pressure,  and  surface 

tension.  It  may  be  expressed: 


u 

m 


f(V  V V 


(4.12) 


The  compressibility  of  the  homogenized  pore  fluid  therefore  can  be 
evaluated  using  the  following  equation. 


Q 


\9u  3e  3u  3e  3t  3 cv) 
a v w v c v/ 


(4.13) 


If  the  form  of  Equation  4.13  is  known,  the  compressibility  of  the 
homogenized  pore  fluid  can  be  evaluated  by  knowing  the  changes  in  air 
pressure,  water  pressure,  and  tension  due  to  unit  change  in  volumetric 
strain.  For  the  case  in  which  the  air  bubbles  are  occluded,  Eq.  4.13 
reduces  to: 
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(4.14) 


3u 

w 


For  this  important  practical  case,  the  compressibility  of  the 
homogenized  pore  fluid  can  be  evaluated  using  Boyle's  and  Henry's  Law. 
The  details  will  be  given  later. 

For  clays  with  low  degrees  of  saturation,  open  channels  are  likely 
to  be  present.  The  behavior  in  this  case  is  extremely  complicated,  and 
no  consideration  has  been  given  to  it  during  this  study. 


ASSUMPTIONS  FOR  DERIVATION  OF  THE  COMPRESSIBILITY  OF  THE 
HOMOGENIZED  PORE  FLUID 

(1)  Boyle's  Law 

Boyle's  Ixiw  states  that  the  product  of  pressure  and  volume 
of  a gas  is  constant  under  constant  temperature  conditions: 

(2)  Henry's  Law 

Henry's  Law  states  that  at  a constant  temperature,  the  weight 
of  gas  which  can  be  dissolved  in  a given  volume  of  liquid  is 
directly  proportional  to  the  gas  pressure. 

Let  V.  be  the  volume  that  would  be  occupied  by  the  dissolved  gas 
d 

if  it  was  extracted  from  the  liquid  and  compressed  to  the  pressure  acting 
on  the  fluid.  According  to  Boyle's  and  Henry's  Law  this  volume  will  be 
constant  and  thus  independent  of  the  pressure.  is  only  dependent  on 

the  volume  of  water  and  can  be  calculated  as  follows: 

V = H V (4.15) 

d w 

where  V is  the  volume  of  water,  H is  the  coefficient  of  solubility.  At 
w 

20°C,  H - 0.02. 

According  to  Dalton's  division  law,  the  saturated  water  vapor 
pressure  in  the  free  air  does  not  obey  Boyle’s  Law.  However,  the  satu- 
rated water  vapor  pressure  is  usually  very  small  (Schuurman, 1966)  and  the 
influence  can  be  disregarded. 
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According  to  Henry’s  Law,  the  volume  of  air  dissolved  in  the  water 
is  proportional  to  the  volume  of  water.  The  rate  at  which  the  air  dis- 
solves in  the  water  depends  on  the  air  pressure.  The  time  t necessary  to 
dissolve  the  air  for  a unit  change  in  pressure  was  shown  by  Beek  (1963) 
to  be 


t 


(4.16) 


where  is  the  diffusion  coefficient  and  r is  the  radius  of  the  bubble. 
At  20°C,  is  equal  to  10  5 cm2/sec.  For  small  air  bubbles,  the  time 
necessary  to  dissolve  the  air  is  very  small.  Therefore  it  can  be  assumed 
that  the  air  dissolves  in  the  water  instantaneously  and  that  the  water  is 
always  saturated  with  dissolved  air. 

As  a result  of  the  surface  tension,  T , the  air  pressure  and  the 

c 

water  pressure  of  an  air-water  mixture  are  not  the  same.  The  value  of 

u and  u are  related  by 
w a 

u = u - u (4.17) 

w a c 


where  u,  is  the  capillary  pressure  due  to  surface  tension.  uc  is  usually 
ejq>ressed  in  terms  of  the  meniscus  radii,  r^,  as  shown  in  Fig.  4.2. 


(4.18) 


In  a mixture  of  air  and  water,  the  bubbles  do  not  necessarily  all  have 
the  same  shape  and  size.  Therefore,  it  is  assumed  that: 


u 

c 


(4.19) 


where  Rc  is  the  average  capillary  radius  for  the  mixture  of  air  and 

water.  Generally,  for  soils  with  high  degrees  of  saturation,  the  effect 

of  u is  small, 
c 
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FIG.  4.2  DOUBLY  CURVED  INTERFACE 
BETWEEN  AIR  AND  WATER 
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DERIVATION  OF  THE  COMPRESSIBILITY  OF  A MIXTURE  OF  AIR  AND  WATER 


If  e is  the  volume  of  free  air  and  e is  the  volume  of  dissolved 
a s 

air  in  the  air-water  mixture,  then  it  follows  from  Boyle's  and  Henry's 
Law  that 


(ea+es)(pa+ua)  = 


(eaO  + es)(pa+UaO) 


(4.20) 


where  p is  atmospheric  pressure,  e n is  the  initial  air  volume  and  u „ 
a au  aO 

is  the  initial  air  gauge  pressure.  Equation  4.20  can  also  be  expressed 

as: 


u = 
a 


(e  - e )p  + (e  + e )u  _ 
aO  a a aO  s aO 

e + e 
a s 


(4.21) 


Using  Equation  4.21  and  Equation  4.17,  the  pore  water  pressure  may 
be  expressed  as: 


u = 
w 


(e  - e )p  + (e  . + e )u  „ 2T 
aO  a a aO  s aO  c 

e + e R 

as  c 


(4.22) 


The  compressibility  of  the  air-water  mixture  is 


1 

Q 


de 
\ 

du 


de 
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(l  + e0)  duw 


(4.23) 


where  eQ  is  the  initial  void  ratio.  From  equation  4.22 


du 
v 

de 


(e  . + e )p  + (e  . + e )u  . 2T  dR 
aO  s *a  aO  s aO  c c 

(e  + e ) 2 R 2 de 

as  c a 


(4.24) 


In  terms  of  degree  of  saturation  S and  void  ratio  e 


du 
w 

de 


^'Vo^Vo*  (ua0+pa) 
(e  - Se  + HSe)  2 


2T  dR  „ 
c c £ 

_ 2 dS  e 


(4.25) 
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Substituting  Eq.  4.23  into  Eq.  4.25,  the  compressibility  of  the  homo- 
genized pore  fluid  can  be  expressed  as: 


1 

Q 


(l  + e0) 


(eQ-SOeO  + HSOeO)(uaO  + Pa) 


(e  - Se  + HSe) 


2T  dR  „ 
c c £ 

D 2 dS  e 


(4.26) 


COMPARISON  BETWEEN  THEORIES  AND  TEST  DATA 

The  relationships  described  in  the  previous  sections  can  be  used  bo 
predict  the  change  in  the  degree  of  saturation  of  a sample  resulting  from 
application  of  a back  pressure. 

If  the  surface  tension  term  in  Eq.  4.26  is  neglected,  the  equation 
reduces  to  a form  similar  to  that  suggested  by  Hilf  (1948)  and  Skempton 
and  Bishop  (1954) . The  comparison  between  test  data  obtained  by 
Mitchell  et  al,  (1965)  and  the  results  predicted  by  Hilf's  and  Skenpiton 
and  Bishop's  equation  is  shown  in  Fig.  4.3. 

Schuurman  (1966)  took  surface  tension  into  account  by  assuning  that 
the  air  bubbles  were  all  of  equal  size  and  spherical  in  shape.  By  this 
assumption,  the  value  of  Rc  in  Eq.  4.19  would  be  equal  to  the  radius  of 
the  air  bubbles.  Taking  the  initial  radius  of  air  bubbles  as  RQ,  the 
radius  after  a change  in  void  ratio  can  be  expressed  as 


Substitute  this  expression  into  Eq.  4.25,  the  compressibility  of 
this  air-water  mixture  can  be  calculated.  A comparison  between  test 
data  obtained  by  Mitchell  et  al.  (1965)  and  the  results  predicted  by 
Schuurman 's  equation  is  shown  in  Fig.  4.3. 

It  seems  likely  that  the  discrepancies  between  experimental  data 
and  the  results  predicted  by  Hilf's  and  Skempton  and  Bishop's  equation 
are  mainly  due  to  the  neglecting  of  the  effect  of  surface  tension.  The 
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Degree  of  Soturotion,  S (%) 


Bock  Pressure,  u^kg/cm2) 


FIG.  4.3  CALCULATED  AND  MEASURED  VARIATIONS  OF 

DEGREE  OF  SATURATION  WITH  BACK  PRESSURE 

(Data  from  Mitchell,  et  al.,  1965) 
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discrepancies  between  the  experimental  data  and  the  results  predicted  by 
Schuurman's  equation  are  probably  due  to  the  assumption  that  all  the 
bubbles  are  spherical  and  uniform,  and  that  the  total  number  of  bubbles 
does  not  change  during  the  compression  of  the  soil. 

A further  consideration  relative  to  surface  tension  effects  is 
that  the  capillary  pressure  does  not  necessarily  arise  only  from  the 
bubbles.  It  can  also  arise  from  the  other  water  surfaces  exposed  to  air 
at  the  specimen  boundaries.  Therefore  it  seems  more  reasonable  to  use 
an  empirically  determined  average  capillary  radius  R^  to  simulate  capil- 
lary effects. 

It  has  been  found  that  the  calculated  results  agree  quite  well  with 
the  test  by  Mitchell  et  al.  (1965)  if  it  is  assumed  that  the  average 
capillary  radius  R^  and  the  degree  of  saturation  are  related  as  follows: 


(4.28) 


where  R is  the  capillary  radius  at  saturation.  The  comparison  is  shown 
in  Fig.  4.4.  is  the  lowest  degree  of  saturation  at  which  the  water 

begins  to  flow  freely.  may  be  found  from  permeability  tests  for 

samples  with  different  degrees  of  saturation.  For  the  clay  tested  by 
Mitchell  et  al.,  was  found  to  be  zero.  It  was  also  found  that  the 
relation  between  the  water  pressure  and  the  degree  of  saturation  was 
approximated  best  by  adopting  a value  of  2TC/RCS  = 5 psi.  At  20°C  the 
value  of  Tc  is  74  x 10”6  kg/cm.  The  value  of  2TC/RCS  = 5 psi  corres- 
ponds to  an  effective  value  of  R^g  equal  to  4.23  x 10~4  cm. 

In  the  case  of  the  New  Melones  Dam  core  material,  sufficient  data 

were  not  available  to  determine  S,  and  2T  /R  . Therefore,  for  purposes 

f c cs 

of  approximating  the  compressibility  of  pore  fluid,  it  was  assumed  that 
both  2Tc/Rcs  and  were  equal  to  zero. 

FLOW  BEHAVIOR  OF  HOMOGENIZED  FLUID 

For  an  isotropic  material,  Darcy's  Law  can  be  ejq>ressed  as 
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(4.29) 


£*■) 


where  z is  the  elevation  and  v is  the  superficial  velocity  of  the  pore 
fluid  relative  to  the  soil  skeleton. 

For  the  case  of  clay  with  occluded  bubbles,  u is  equal  to  u , and 

m w 

Darcy's  Law  may  be  expressed  as 


v = k V 


(€H 


(4.30) 


For  the  case  of  clay  with  a low  degree  of  saturation,  it  is  not 
known  whether  Eq.  4.30  is  valid.  However,  such  conditions  are  not  likely 
to  arise  for  compacted  clay  cores  in  embankment  dams. 
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CHAPTER  5 

FINITE  ELEMENT  ANALYSIS  OF  CONSOLIDATION  OF  SATURATED  CLAY 
INTRODUCTION 

The  equations  governing  the  three-dimensional  consolidation  of 
elastic  media  were  developed  by  Biot  (1941,  1955,  1956) . However,  until 
recently  Biot's  Theory  had  only  been  applied  to  problems  with  simple 
geometry  and  material  properties.  The  solutions  of  one-dimensional  con- 
solidation problems  were  achieved  by  expansion  in  series,  while  integral 
transform  techniques  were  used  in  the  solution  of  two-dimensional 
problems  by  Biot  (1941),  Gibson  et  al.  (1970),  McNamee  and  Gibson  (1960a, 
b)  and  Booker  (1973) . 

The  development  of  the  finite  element  method  provides  a necessary 
tool  for  the  treatment  of  general  consolidation  problems,  with  no  restric- 
tions on  geometry,  material  properties  or  boundary  conditions.  Sandhu 
and  Wilson  (1969)  used  Gurtin's  variational  method  (1964)  to  convert  the 
field  equations  to  a set  of  integral  equations  and  to  formulate  a func- 
tional for  finite  element  analyses.  Later  Hwang  and  Morgenstern  (1971) 
presented  a finite  element  formulation  based  on  the  method  of  weighted 
residuals.  Yokoo  et  al.  (1971a,  b,  c)  investigated  the  use  of  discon- 
tinuous functions  to  improve  the  accuracy  of  the  calculation  of  initial 
pore  pressures. 

During  the  course  of  this  study,  Biot's  equations  were  reviewed, 
and  finite  element  formulation  of  these  equations  were  developed  using 
Galerkin's  process.  An  incremental  finite  element  formulation  was  also 
developed  for  analysis  of  materials  with  non-linear  stress-strain 
behavior.  Several  consolidation  problems  involving  linear  elastic  satu- 
rated materials  were  analyzed  and  the  results  were  compared  with  analy- 
tic solutions. 
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FIELD  EQUATIONS 


The  equations  governing  the  consolidation  of  a saturated  soil  with 
an  elastic  skeleton  are 

(a)  Equilibrium  equation.  Only  changes  in  stresses  and  body 

forces  are  considered.  The  geostatic  body  forces  are  con- 
sidered as  initial  stresses 


[3] {a}  = {fb> 


(5.1) 
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(b)  Continuity  equation.  For  a completely  incompressible  pore 
fluid,  the  amount  of  volumetric  strain  of  soil  equals  the 
amount  of  water  expelled  from  the  soil  skeleton.  For  a 
compressible  pore  fluid,  the  volume  decrease  of  the  soil  is 
the  sum  of  the  volume  expelled,  plus  the  decrease  in  volume  of 
the  fluid  within.  The  continuity  equation  can  therefore  be 
written  as 


W)T  (v)  ♦ (V>T  ^ . i |H  . o 


where  {V} 


/ J.  _i  Jl\ ' 

\ 3x  3y  3z  / 
,T 


(5.4) 

(5.5) 


{v}  * (v^  Vy  vz) * is  the  superficial  velocity  of  the  pore 
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T 

fluid  relative  to  the  soil  skeleton,  {w}  = (w  , w , w ) is 

x y z 

the  excess  pore  pressure,  and  1/Q  is  the  compressibility  of 
pore  fluid. 

(c)  Total  stresses  {a}  are  the  sum  of  effective  stresses  {0}'  and 
pore  pressure  u. 

{a}  = {a}'  + u{a}  (5.6) 

where  {a}  = (1,  1,  1,  0,  0,  0)T  (5.7) 

(d)  The  effective  stress-strain  relationship  can  be  expressed  as 

{o}'  = [CE]{£>  (5.8) 

T 

where  {e}  = (E  E E Y Y Y ) is  the  vector  of 
x y z yz  zx  xy 

strain  components  and  the  strains  {e}  are  defined  as  follows: 


{e}  = mT  {w} 


(5.9) 


where  {w}  is  displacement  vector,  [3]T  is  the  transpose  of 

[9]  defined  by  Eq.5.3.  [cE]  is  the  matrix  of  elastic  constants 

relating  effective  stress  and  strain.  For  an  isotropic 
E 

elastic  body,  [C  ] has  the  form 
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(5.10) 


where  A is  one  of  Lame’s  constants,  and  G is  the  shear  modulus, 
which  may  be  related  to  the  drained  values  of  Young's  modulus 
E'  and  Poisson’s  ratio  V'  as  follows 
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\ = 


E' V' 
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(f)  Darcy's  Law.  For  many  soils,  it  can  be  assumed  that  Darcy's 
Law  is  valid.  The  change  in  superficial  velocity  {v}  rela- 
tive to  the  soil  skeleton  can  be  related  to  the  excess  pore 
pressure  by  the  equation: 


{v}  = — [k]  {V}  • u 


(5.11) 
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where  {v}  and  {V}  are  as  previously  defined,  y^  is  the  unit 
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Substituting  Eqs.  5 

.6,  5 

.8,  5.9 

into  Eq.  5.1, 

the  equilibrium 

equation  can  be  written  as  follows: 

m[cEimT{ 

w}  + [3]{a}u  = 

<fb> 

(5.13) 

Substituting  Eq.  5.11  into  Eq.  5.4,  the  continuity  equation  becomes 


- {V}T [k] {7}u  + {V}T  M - | * 0 (5.14) 

To  completely  define  the  problem  both  elastic  and  flow  boundary  conditions 
must  be  specified.  For  the  elastic  boundary  conditions,  it  is  assumed 
that  part  of  the  boundary  surface  (ST)  is  subjected  to  known  applied 
tractions,  while  the  remainder  of  the  surface  (S^)  is  subjected  to  speci- 
fied displacements.  For  the  flow  boundary  condition,  it  is  assumed  that 
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part  of  the  boundary  surface  (Sq)  is  subjected  to  a specified  flow 
velocity,  while  the  remainder  of  the  surface  (S^)  is  subjected  to  known 
pore  pressures.  The  boundary  conditions  may  be  written  in  the  form 


M = {0}  on  SD  for  t ^ 0 

[N]  {a}  = {Tf}  on  ST  for  t ^ 0 

u = u on  S for  t > 0 

P 

{v}T{n}  = q on  S for  t > 0 

q 


(5.15) 


T 

where  {n}  = (n^  n^  nz)  is  the  outward  normal  to  the  boundary  surface,  [N] 

is  defined  as  follows 
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(5.15a) 


and  (~)  indicates  a prescribed  quantity. 

To  completely  describe  the  problem,  the  initial  conditions  must 
also  be  defined.  When  loads  are  applied,  the  pore  fluid  is  expelled  at 
a finite  rate.  Thus  the  decrease  in  volume  of  a soil  element  must  be 
equal  to  the  decrease  in  volume  of  the  pore  water,  so  that 


e 

v 


m 1 ^u 

= q at 


at  t 


(5.16) 


T 

where  = {a}  {e}  is  the  volumetric  strain,  and  {a}  and  {e}  are  as 
previously  defined. 


For  an  elastic  soil  skeleton,  the  undrained  solution  can  also  be 
obtained  by  using  a bulk  modulus  equal  to  the  sum  of  the  drained  bulk 
modulus  of  the  soil  skeleton  and  the  bulk  modulus  of  the  pore  fluid, 
and  a shear  modulus  equal  to  the  drained  shear  modulus  of  the  skeleton. 
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FINITE  ELEMENT  DISCRETIZATION 


In  Biot's  theory  of  consolidation,  the  field  variables  {w}  and  u 
within  an  element  can  be  discretized  spatially  as  follows: 


w = (N  ) {w  } 
x w x 

w = (N  ) {w  }€ 
y w y 

w = (N  ) {w  }€ 
z w z 

u = (Nu) {u}e 


(5.17) 

(5.18) 

(5.19) 

(5.20) 


where  {w  }e,  {w  }e,  {w  }6  and  {u}6  are  the  values  of  displacements  and 
x y z 

pore  pressure  at  the  nodes  of  this  element,  (N^)  is  the  shape  function 
for  displacement,  and  (N  ) is  the  shape  function  for  pore  pressure. 


Eqs.  5.17,  5.18,  5.19  can  be  combined  as: 


{w}  = [<J>w]{w} 


(5.21) 


where 
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By  using  Galerkin's  process,  Eqs.  5.13  and  5.14  can  be  approximated 
as  a set  of  integral  equations. 
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n-iyv  n“l  jy 

C i4>wHV  « + l / tv  {V  dv 

Js  •'v 


dV 


(5.24) 
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f ^-(N  )T{V}T[k]{V}  u dv  + £ f (N  )T{V}T  -2-  M dV 
n=l  •' v ru)  n=i  Jv 

E / 5,n / K * • / <*.'*  5 - 

n=-l  •'v  * yc 


where  N is  the  total  number  of  elements.  The  integrals  are  evaluated 
over  each  element  and  summed  for  the  entire  body.  Substituting  Eqs. 

5.17  - 20  into  Eqs.  5.24  and  5.25,  the  field  variables  u and  {w}  for  an 
element  can  be  expressed  in  terms  of  the  nodal  values  {u}  and  {w}  and 
Eqs.  5.24,  5.25  yield  a set  of  simultaneous  equations  (5.26,  5.27)  with 
a finite  number  of  nodal  variables  {u}n  and  {w}n,  where  the  superscript 
n indicates  the  vector  of  nodal  values. 

[K] {w}n  + [L]{u}n  = {F}  + {Ffa}  (5.26 


- [H]  {u}n  + [G]  ^ {w}n  - (E]  ± {u}n  = (D 


N r 

where  [K]  = £ / 

n=l  Jv 


(B]T[CE](B)  dV 


IB]  = [33  t4>  1 

W 


N r 

E / <*» 

n=l  Jv 


] [31 (a> (N  ) dv 
u 


[G1  = 2^  j (Mu)T{V}T[^wl  dV  - [L]T 


t»]  - £ f , )T{V}'iT1c1  {V}  (N  ) dV 
n-1  Jv  Tu>  u u 

,E1  * J(  / 5 ‘“/"V  * 
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(5.34) 


{F} 


{F  } = 

D 


{r} 


l 
s/w- 
/ 


) {Tf } dS 


] { f . } dV 
D 


(N  ) q dS 
u 


(5.35) 


(5.36) 


NUMERICAL  INTEGRATION 

Eq.  5.27  contains  terms  ^ {u}n  and  ^ {w}n.  Therefore  an  integra- 
tion over  time  is  necessary.  It  is  assumed  that,  to  a sufficient  degree 
of  approximation: 


i: 


{U}n  dt  « 0 At{u}t  + (1  -e)&t  {u)t 


(5.37) 


and 


/ 


{w}n  dt  = 0 At{w}”  + (1-6)  At  {w}" 

1 0 


(5.38) 


where  0 is  a factor  whose  value  may  vary  from  0 to  1.  A value  of  0 = -? 

is  equivalent  to  the  assumption  of  a linear  variation  in  time  and  thus 
the  central  difference  approximation  (the  Crank-Nicholson  Method) . The 
approximation  6=0  leads  to  the  fully  explicit  method,  and  the  approxi- 
mation when  0 = 1 leads  to  the  fully  implicit  method.  When  0 < 1, 
it  can  be  shown  that  numerical  integration  is  unconditionally  stable 
(Booker  and  Small,  1974). 

Integrating  Eq.  5.27  with  the  aid  of  Eqs.  5.37,  5.38  leads  to 

-At[H]  ( (1  - 0)  {u}"  + 0{u}"  )+  IGj({w}"  - M"  )- 

v 'o  tl/  1 'o' 

[E]  ({u}"  - {u}"  )-  At  ("(1-6)  {r}"  +0  {I*}"  ) (5.39) 

V 'o/  \ 'o  1' 
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or 


[G]  {w}”  - ([H]0  At  + [E]  ) {u}”  = [G]{w}"  + 

1 V 7 1 t0 

(iHHi-ejAt  - ie]  ) {u}"  + At  f (l- ejtn"  + e {r}J  ) (5.40) 

7 0 v Si7 

Combining  Eq.  5.40  and  Eq.  5.26  and  observing  that 
[G]  = [iijT 

it  is  found  that 

tK]  [L1  1 = (F  + M 

[L]t  -[H]6  At-  [E]J  ( ujt!  ( 0Atr  j t 

l 

I (1  - 0) At  {r}  + ( [H] (1  - 0)At  - [E] ) {u}  + [L]T{w> 

0 0 0 

Eq.  5.41  shows  that  if  the  solution  is  known  at  tQ,  then  it  can  be  « 

found  at  t^,  and  the  solution  can  thus  be  obtained  by  marching  forward. 

Formulation  for  an  Incremental  Effective  Stress-Strain  Law 

It  is  often  possible  to  describe  the  behavior  of  an  inelastic 
material  by  an  incremental  stress-strain  law  of  the  form 

{a}  = [Ct]{e}  (5.42) 

where  {o}  = ~ and  {e}  * are  the  stress  and  strain  rates,  respectively. 

lCt]  is  a matrix  of  tangent  values  which  are  functions  of  the 
current  state  of  stress.  In  this  study,  am  elasto-plastic  hardening 
material  with  an  associated  flow  rule  was  used.  The  matrix  [Ctl  is  de- 
rived in  Appendix  A. 
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| (5.41) 


For  such  a material  it  can  be  shown  that  Eq.  5.13  must  be  re- 
placed by 

[3]  [C  ] [9 ]T{w}  + [31 {a}  u=  {f.  } (5.43) 

t b 

Eq.  5.14  remains  valid  for  an  increment  of  displacement  and  pore 
pressure,  however. 

The  derivation  of  the  finite  element  formulation  is  quite  similar 
to  that  described  previously,  and  leads  to  the  equations 

[L]  I l Aw 

- [H]0  At  - [E]J  ( u. 

fcl 

AF  + AFb  \ j [L]{u}t 

x+  d-e)At  (r}t0)  ( ( [hi  (i  - e>  At  - [ei){uJ 

0 

where  [ Kfcl  = J [B ] T (C^ ] [Bl  dV  is  the  average  value  of  the 'incremental 

*v 

tangent  stiffness  matrix  for  this  increment. 

Eq.  5.44  is  a nonlinear  equation.  An  iteration  process  is  neces- 
sary in  order  to  obtain  an  accurate  solution. 

Iteration  Scheme 

Embankment  construction  is  simulated  using  the  technique  developed 
by  Clough  and  Woodward  (1967) . For  each  increment  of  loading,  the  stiff- 
ness matrix  is  recalculated  using  the  current  material  properties  and 
the  current  state  of  stress.  For  each  increment  of  loading,  an  initial 
strain  iteration  scheme  is  used  to  solve  the  approximating  equation 
(5.44)  . 

EXAMPLES 

To  demonstrate  the  accuracy  and  applicability  of  the  finite  element 


| (5.44) 


[Kt] 

[LlT 


6At  (rlt 
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method  presented  in  this  study,  a number  of  problems  were  analyzed  and 
the  results  were  compared  with  analytic  solutions. 

(1)  One-dimensional  consolidation  of  a homogeneous  layer.  Fig. 

5.1  shows  the  isochrones  in  a single  layer  with  an  impervious 
base.  It  can  be  seen  that  the  finite  element  solution  approxi- 
mates the  exact  solution  of  this  problem  very  closely. 

(2)  One-dimensional  consolidation  of  contiguous  layers.  The 
second  problem  analyzed  was  that  of  two  contiguous  layers 
with  the  same  compressibility,  the  upper  one  being  four  times 
as  permeable  as  the  lower  one.  The  results  obtained  by  five 
analyses  are  shown  in  Fig.  5.2.  They  are: 

1.  Schmidt's  graphical  procedure  (1924) 

2.  Euscher's  analog  computer  solution  (1965) 

3.  Harr's  finite  difference  solution  (1967) 

4.  Christian's  finite  element  solution  (1970) 

5.  The  finite  element  developed  solution  developed 
during  this  study. 

It  can  be  seen  that  there  is  considerable  spread  among  these 
results.  Because  of  the  approximations  it  involved,  Luscher 
expected  his  result  to  be  too  slow.  Harr's  explicit  finite 
difference  procedure  does  not  consider  the  compressibility  of 
the  nearby  material  at  the  interface  between  two  materials, 
and  it  therefore  leads  to  too  rapid  consolidation.  Christian's 
finite  element  solution,  Schmidt's  graphical  procedure,  and 
the  finite  element  solution  developed  during  this  study  are  in 
good  agreement,  and  are  believed  to  be  accurate.  However,  it 
should  be  emphasized  that  all  the  solutions  to  this  problem 
are  approximate. 

(3)  A strip  footing  on  a layer  of  finite  thickness.  The  finite 
element  mesh  used  for  the  analysis  is  shown  in  Fig.  5.3.  The 
half  width  of  footing  is  a.  The  soil  was  assumed  to  be  linear 
and  elastic,  with  a shear  modulus,  G,  and  a Poisson's  ratio,  v' . 
The  finite  element  solution  is  compared  with  the  analytic 
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Depth  Ratio,  Z/H 


Excess  Pore  Pressure,  u/u0 


FIG.  5.2  CALCULATED  PORE  PRESSURE  ISOCHRONES  FOR 
CONTIGUOUS  LAYERS  WITH  DIFFERENT 
PERMEABILITIES 
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4a 


FIG.  5.3  FINITE  ELEMENT  MESH  FOR  A STRIP 
FOOTING  ON  A FINITE  LAYER 
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solution  developed  by  Gibson,  et  al.  (1970) . In  Fig.  5.4, 
the  displacement  is  plotted  against  time  factor.  It  can  be 
seen  that  the  results  of  the  finite  element  analysis  are  in 
good  agreement  with  the  closed  form  solution. 

(4)  One-dimensional  incremental  deposition  to  simulate  the  time- 
dependent  deposition  and  simultaneous  consolidation  of  a 
layer  of  material;  elements  are  added  to  the  mesh  at  inter- 
vals of  time.  As  shown  in  Fig.  5.5  element  no.  2 is  placed 
on  top  of  element  no.  1 and  element  no.  3 is  placed  on  top  of 
element  no.  2 to  simulate  the  increase  of  the  layer  thickness 
with  time.  The  rate  of  deposition  is  m (units  of  length/time). 
The  time  factor  used  for  this  problem  is  T = m2t/cv.  When  T 
is  infinite,  deposition  is  instantaneous,  and  no  consolidation 
occurs  during  deposition.  When  T is  small,  it  represents 
slow  deposition,  and  appreciable  consolidation  occurs  during 
deposition.  Figs.  5.6  and  5.7  show  the  pore  pressure  dis- 
tributions calculated  for  different  values  of  the  time  factor, 
T,  for  layers  deposited  on  both  permeable  and  impermeable 
bases.  The  finite  element  solution  shows  good  agreement  with 
Gibson's  (1958)  analytic  solution  for  both  cases. 
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Dimensionless  Displacement, 


rime  Factor,  T * cvt/H2 

I0*4  I O'3  IO"2  ICT1  I 10 


_FIG.  5.4  SURFACE  DISPLACEMENT  OF  A STRIP  FOOTING  ON 
A LAYER  OF  FINITE  THICKNESS 
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Pore  Pressure  Rotio,  u/yH 


FIG.  5.6  PORE  DISTRIBUTIONS  FOR  DEPOSITION  ON  AN 
IMPERMEABLE  BASE 
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Height  Rotio,  Z/H 


CHAPTER  6 


FINITE  ELEMENT  ANALYSIS  OF  CONSOLIDATION  OF  PARTLY  SATURATED  CLAY 


) 


1 


INTRODUCTION 

As  discussed  in  Chapter  4,  a three-phase  unsaturated  soil  can  be 
treated  phenomenologically  as  a two-phase  saturated  material  with  a com- 
pressible pore  fluid.  This  means  that  a partly  saturated  soil  can  be 
analyzed  using  the  theory  of  consolidation  developed  in  the  pervious 
chapter,  provided  the  correct  values  of  permeability  and  compressibility 
are  used.  The  permeability  and  compressibility  of  partly  saturated  clay 
were  discussed  in  Chapters  3 and  4 . 

In  this  chapter,  the  finite  element  formulation  is  derived  for  an 
unsaturated  soil  with  a pore  fluid  which  contains  occluded  bubbles. 

This  important  practical  case  frequently  occurs  where  the  core  of  an  earth 
dam  is  compacted  to  a water  content  in  the  neighborhood  of  optimum.  A 
number  of  cases  were  analyzed  to  demonstrate  the  applicability  of  this 
theory . 

FIELD  EQUATION 

The  equations  governing  the  consolidation  of  a saturated  material 
with  a homogenized  compressible  pore  fluid  are: 

(1)  the  equilibrium  equation, 

(2)  the  effective  stress  equation, 

(3)  the  continuity  equation, 

(4)  Darcy's  Law,  and 

(5)  the  effective  stress-strain  relationship 

These  equations  have  been  derived  in  Chapter  5 for  a material  satu- 
rated with  water.  Partly  saturated  materials  may  be  considered  to  be 
saturated  with  a homogenized  compressible  pore  fluid,  and  the  equations 
governing  the  consolidation  behavior  are  the  same  as  those  developed  in 
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Chapter  5,  except  that  the  pore  water  pressure  is  now  replaced  by  the 
homogenized  pore  fluid  pressure. 

As  discussed  in  Chapter  4 the  homogenized  pore  fluid  pressure  is 
equal  to  the  pore  water  pressure  for  clays  compacted  wet  of  optimum. 
Therefore,  in  this  case,  the  equations  governing  the  consolidation 
behavior  are  precisely  those  developed  in  Chapter  5. 

FINITE  ELEMENT  FORMULATION 

The  derivation  of  finite  element  formulation  is  the  same  as  that 
given  in  Chapter  5 and  will  not  be  repeated  here. 

The  matrix  is: 


Af  + AF. 


[L]{u}t 


0At{r}  + (i  - 8)At{r}t  j i ([h] (l- 0) At  - [e]){u} 


(H]  = £ f 7^(N  )T{V}' 
n-1  Jv  ru> 

■ E / | <»u>t,nu 

n*l  J v 


)T{V}T[k]{V} (Nu)  dv 


(6.9) 


(6.10) 


In  Eq.  6.9  and  Eq.  6.10,  the  permeability  and  compressibility  are 
functions  of  both  degree  of  saturation  and  void  ratio.  As  described  in 
Chapter  3 and  Chapter  4,  Eq.  3.2  and  Eq.  4.26  can  be  employed  to  calcu- 
late the  appropriate  values  of  permeability  and  compressibility. 

The  computer  program  developed  based  on  this  theory  is  described 
in  Appendix  E.  Isoparamatric  elements  with  four  to  eight  nodes  (Bathe 
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and  Wilson,  1973)  are  used  in  the  computer  program.  The  advantage  of 
using  this  type  of  element  for  analyses  of  consolidation  are: 

(1)  The  numerical  difficulties  caused  by  sudden  variations  in 
pore  pressure  gradients  are  reduced  considerably  by  using 
four  point  Gaussian  integration.  This  method  utilizes  higher 
order  elements  capable  of  better  approximating  variations  in 
pore  pressure  gradients. 

(2)  Any  number  of  nodes  from  four  to  eight  can  be  specified  for 
each  element.  Triangular  shaped  elements  can  be  formed  by 
superimposing  nodes.  This  allows  considerable  flexibility 
in  mesh  layout  for  the  complicated  geometry  of  a zoned  dam. 

EXAMPLES 

Two  examples  are  given  in  this  chapter  to  demonstrate  the  applica- 
tion of  this  theory,  to  check  the  accuracy  of  numerical  calculation,  and 
to  compare  the  consolidation  behavior  of  saturated  and  partly  saturated 
soil. 

As  a first  example,  the  one-dimensional  consolidation  test  for  a 
linear  elastic  material  with  a constant  permeability  and  a constant  pore 
fluid  compressibility  was  considered. 

This  problem  gives  some  insight  into  the  consolidation  behavior  of 
unsaturated  soils  and  has  the  advantage  that  it  was  possible  to  deduce 
an  analytic  solution,  so  that  the  accuracy  of  the  finite  element  approxi- 
mation can  be  checked. 

When  the  soil  is  saturated  with  water,  the  compressibility  of  pore 
fluid  (1/Q)  is  very  small,  and  the  equations  reduce  to  the  familiar 
one-dimensional  consolidation  equation.  The  solution  is  the  well  known 
Terzaghi's  solution.  When  the  soil  is  con^letely  dry,  the  value  of  pore 
fluid  compressibility  is  infinite,  and  no  consolidation  takes  place.  For 
partly  saturated  soils,  the  value  of  Q lies  between  zero  and  infinity, 
and  some  behavior  intermediate  between  the  extremes  discussed  above  may 
be  expected. 
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For  the  purpose  of  comparing  the  consolidation  behavior  of  soils 
with  different  degrees  of  compressibility  of  pore  fluid  it  is  convenient 
to  define  a dimensionless  parameter  n : 


(6.11) 


where  D'  is  the  constrained  modulus  of  the  soil  skeleton. 

The  one-dimensional  consolidation  behavior  was  calculated  with  a 
range  of  T)  from  1 to  20000  using  the  finite  element  method  described 
previously.  The  time-settlement  curves  are  shown  in  Fig.  6.1,  and  the 
pore  pressure  dissipation  curves  are  shown  in  Fig.  6.2.  It  may  be  seen 
that  the  results  obtained  from  the  finite  element  analyses  agree  well 
with  those  obtained  from  a theoretical  solution  given  in  Appendix  B. 

The  results  indicate  that  a significant  initial  settlement  occurs 
for  materials  with  highly  compressible  pore  fluid.  As  time  goes  on,  the 
rate  of  settlement  is  retarded  by  the  effect  of  the  compressibility  of 
the  pore  fluid.  These  characteristics  are  often  observed  in  the  field. 
The  same  tendency  has  been  found  in  10-in. -dia.  oedometer  tests  on  undis- 
turbed sanqples  of  compacted  clay  fill  from  Derwent  Dam  (Barden,  1965) . 

The  material  properties  used  in  the  above  example  are  an  over- 
simplification of  the  properties  of  a real  soil.  For  the  second  example, 
the  one-dimensional  consolidation  behavior  was  studied  using  the  proper- 
ties of  the  soft  core  material  from  New  Melones  Dam. 

It  was  assumed  that  the  soil  is  an  elastic-plastic  material  with 
stress-strain  relationships  of  the  type  described  in  Chapter  2.  The 
summary  of  parameters  for  the  stiff  core  material  of  New  Melones  Dam  is 
given  in  Table  6.1.  The  determination  of  these  parameters  is  discussed 
in  Appendix  D. 

The  permeability  of  the  pore  fluid  and  the  permeability  of  the  soil 
vary  with  degree  of  saturation  and  void  ratio,  as  discussed  in  Chapter  3 
and  Chapter  4. 


FIG.  6.1  VARIATION  OF  DEGREE  OF  CONSOLIDATION  WITH  TIME  FOR  ONE- 
DIMENSIONAL CONSOLIDATION  OF  AN  IDEALIZED  ELASTIC  MATERIAL 


FIG.  6.2  PORE  PRESSURE  DISSIPATION  WITH  TIME  FOR  ONE-DIMENSIONAL 
CONSOLIDATION  OF  AN  IDEALIZED  ELASTIC  MATERIAL 


The  relationship  between  the  degree  of  saturation,  the  void  ratio 
and  the  permeability  used  for  this  example  is 


k = k S' 
s 


1 + e 


/— 

/ 1 + € 


(6.12) 


where  eQ  is  the  initial  void  ratio. 

Eq.  4.26  was  used  in  this  example  to  calculate  the  compressibility 
of  the  pore  fluid,  and  it  was  assumed  that  the  surface  tension  term  in 
Eq.  4.26  was  negligible.  Thus  the  compressibility  of  the  pore  fluid  was 
expressed  as 


1 

Q 


1 + e„ 


(eo  + HSeo  - Seo> 

(e  + HSe  - Se)  2 


(ua0  + Pa) 


(6.13) 


Two  conditions  were  considered.  The  first  was  a one-dimensional 
consolidation  test  of  a saturated  sample,  and  the  second  was  a one- 
dimensional consolidation  test  on  a partly  saturated  sample  compacted  to 
a degree  of  saturation  of  85%.  The  time-settlement  curves  of  these  two 
samples  were  calculated  in  three  load  increments.  The  results  are  shown 
in  Fig.  6.3. 

During  the  first  load  increment,  from  1 tsf  to  4 tsf,  the  compres- 
sibility of  the  pore  fluid  was  high  in  the  partly  saturated  sample,  and 
the  initial  settlement  was  large  compared  to  the  subsequent  consolidation 
settlement.  During  the  load  increment  from  4 tsf  to  16  tsf,  the  degree 
of  saturation  of  the  partly  saturated  sample  increased,  so  that  less 
difference  between  the  shapes  of  the  consolidation  curves  for  the  satu- 
rated and  partly  saturated  samples  was  observed. 

During  the  load  increment  from  16  tsf  to  40  tsf,  the  unsaturated 
sample  was  sufficiently  compressed  to  become  saturated.  Therefore,  the 
time  settlement  curve  is  essentially  identical  to  that  of  the  saturated 
sample. 
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FIG.  6.3  TIME-SETTLEMENT  CURVE  FOR  THE  NEW  MELONES  DAM  CORE  MATERIAL 


It  is  interesting  that,  for  this  soil,  the  calculations  indicate 
that  the  saturated  sample  would  undergo  more  settlement  than  the  unsatu- 
rated sample  over  a long  period  of  time.  The  reason  for  this  is  that 
although  the  unsaturated  sample  has  a greater  initial  settlement,  con- 
solidation proceeds  at  a slower  rate  because  of  the  compressibility  of 
the  pore  fluid,  and  eventually  the  settlement  of  the  saturated  sample 
exceeds  that  of  the  unsaturated  sample . 
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CHAPTER  7 
NEW  MELONES  DAM 


INTRODUCTION 

The  procedures  for  analysis  of  the  consolidation  of  partly  satu- 
rated soils  which  was  described  in  the  previous  chapters  was  used  to 
analyze  the  pore  pressures  and  the  movements  in  the  core  of  New  Melones 
Dam  during  and  after  construction. 

A number  of  dams  have  been  analyzed  to  determine  probable  movements 
during  construction  (Kulhawy,  et  al.,  1969)  and  to  calculate  the  stresses 
in  the  long  term  condition  after  all  excess  pore  pressures  have  dissi- 
pated (Duncan,  1972) . The  analyses  of  long-term  stresses  have  been  per- 
formed assuming  that  the  dams  were  1t  very  slowly,  and  that  there  were 

no  excess  pore  pressures  at  any  ,f  construction.  Although  Koppula 

and  Morgenstern  (1972)  developed  sin^e  procedures  for  approximate 
analyses  of  the  consolidation  of  dams , it  has  not  been  possible  until  now 
to  analyze  the  stresses  and  movements  in  dams  during  consolidation, 
taking  into  account  nonlinear  stress-strain  behavior  and  variations  in 
permeability  and  compressibility  of  the  pore  fluid  as  the  degree  of 
saturation  increases. 

Hie  analysis  of  New  Melones  Dam  described  in  this  chapter  thus 
provides  insight  into  some  important  questions  relating  to  the  behavior 
of  zoned  dams  during  consolidation  of  the  core.  These  include: 

(1)  What  is  the  nature  of  the  expected  movements  in  a zoned  dam 
during  consolidation  of  the  core? 

(2)  How  do  the  stresses  in  the  dam  change  during  consolidation? 

(3)  Are  the  long-term  stresses,  calculated  assuming  slow  construc- 
tion and  no  excess  pore  pressure,  the  same  as  these  calculated 
taking  account  of  consolidation? 


81 


DESCRIPTION  OF  NEW  ME LONE S DAM 


New  Melones  Dam  is  now  (1976)  being  constructed  by  the  U.  S.  Army 
Corps  of  Engineers  to  create  a multi-purpose  reservoir  on  the  Stanislaus 
River.  It  is  located  about  40  miles  east  of  Stockton,  California,  as 
shown  in  Fig.  7.1.  The  dam  will  have  a maximum  height  of  625  feet  above 
the  stream-bed,  and  a crest  length  of  about  1600  feet.  It  will  be  a 
rolled  rockfill  embankment  with  an  impervious  core,  and  will  have  a total 
volume  of  about  16  million  cubic  yards.  The  gross  storage  capacity  of 
the  reservoir  created  by  the  dam  will  be  about  2.4  million  acre-feet. 

A plan  view  of  the  damsite  and  a longitudinal  section  along  the 
dam  axis  are  shown  in  Fig.  7.2.  The  valley  walls  have  average  slopes  of 
1.2  to  1 (horizontal  to  vertical).  The  abutments  and  stream  bottom  have 
been  stripped  of  soil  overburden,  so  that  the  dam  will  be  founded  on 
relatively  hard  metavolcanic  rocks. 

The  maximum  data  section  is  shown  in  Fig.  7.3.  The  dam  will  be 
built  in  three  increments  - a construction  dike  to  elevation  630  feet, 
an  upstream  cofferdam  to  elevation  808  feet,  and  finally  the  main  dam 
to  elevation  1135  feet.  It  is  expected  that  the  dam  will  be  built  within 
about  three  and  a half  years,  and  that  it  will  be  completed  by  the  end  of 
1978. 

ANALYSIS  PROCEDURE 

The  elastic-plastic  stress-strain  theory  described  in  Chapter  2 
was  employed  to  model  the  effective  stress-strain  relationship  of  the 
compacted  clay  in  the  core  of  New  Melones  Dam.  The  homogenized  pore 
fluid  concept  described  in  Chapter  4 was  used  to  treat  the  three-phase 
partly  saturated  compacted  clay  as  a two-phase  material  with  a compres- 
sible pore  fluid.  The  effective  stress  equation  used  in  this  analysis 
is  thus 

O'  - O - u (7.1) 

ID 

in  which  u^  is  the  homogenized  pore  pressure. 
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FIG.  7.1  PROJECT  LOCATION  MAP 
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FIG.  7.3  MAXIMUM  DAM  SECTION 


As  described  in  Chapter  4,  it  is  reasonable  to  write  Eq.  7.1  as 
Eq.  7.2  for  a clay  core  compacted  to  a high  degree  of  saturation 

a'  = a - u (7.2) 

w 

and  this  equation  was  used  for  the  analysis  of  New  Melones  Dam. 

The  compressibility  of  the  homogenized  pore  fluid  is  a function  of 
void  ratio  and  degree  of  saturation  as  described  in  Chapter  4.  This 
compressibility  can  be  calculated  using  the  following  equation: 


1 

Q 


1 + e. 


<«Q  - Se0  * "Vo1 
(e  - Se  + Hse) 2 


(ua0  + Pa) 


(7.3) 


As  explained  in  Chapter  4,  the  effect  of  surface  tension  was  neglected  in 
this  analysis. 


The  permeability  of  partly  saturated  clays  varies  with  degree  of 
saturation  and  void  ratio,  as  described  in  Chapter  3.  Because  there  were 
not  sufficient  data  for  the  determination  of  the  functions  which  relate 
permeability,  degree  of  saturation  and  void  ratio,  Poisseuille's  equation 
was  used  to  approximate  this  relationship.  According  to  Poiseuille's 
equation,  the  permeability  at  any  degree  of  saturation  is  related  to  the 
saturated  permeability  as  indicated  by  Eq.  7.4. 


k * k S3 
s 
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(7.4) 


PROPERTIES  OF  THE  MATERIALS  IN  NEW  MELONES  DAM 

New  Melones  Dam  contains  a number  of  zones  of  different  materials, 
as  shown  in  Fig.  7.3.  For  purposes  of  the  analysis  described  in  this 
chapter,  the  embankment  cross-section  was  simplified  as  shown  in  Fig.  7.4. 
This  simplified  cross-section  contains  two  materials — the  core  and  the 
shell. 
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Because  variations  in  water  content  and  dry  density  may  occur  dur- 
ing construction  of  the  dam,  two  different  compaction  conditions  were 
considered  for  the  core  material:  a "stiff"  core  (corresponding  to  95% 

relative  compaction  as  determined  by  the  Standard  AASHO  compaction  test, 
and  1%  dry  of  optimum  moisture  water  content) , and  a "soft"  core  (corres- 
ponding to  90%  relative  compaction  and  optimum  moisture  content  at  the 
time  of  placement) . These  are  the  same  conditions  as  considered  by 
Quigley,  et  al.  (1976)  in  their  analysis  of  stresses  and  movements  in 
New  Melones  Dam,  which  was  made  using  hyperbolic  parameters. 

The  procedures  used  to  determine  the  values  of  the  parameters  of 
the  shell  and  the  core  material  for  the  Cam  Clay  Model  are  described  in 
Appendix  D. 

ANALYSIS 

It  was  not  feasible  to  perform  a three-dimensional  analysis  of  the 
New  Melones  Dam.  It  was  assumed  that  a reasonable  approximation  to  the 
behavior  at  the  maximum  section  of  the  dam  could  be  obtained  by  perform- 
ing a plane  strain  analysis  of  a typical  section  at  the  center  of  the 
valley. 

The  mesh  used  for  this  analysis  is  shown  in  Fig.  7.4,  and  the 
computer  program  used  for  the  analysis  is  given  in  Appendix  E. 

The  computer  program  models  incremental  loading,  so  that  the  con- 
struction of  a dam  can  be  simulated.  The  program  also  has  the  capability 
of  simulating  the  filling  of  a reservoir  behind  an  embankment  by  specify- 
ing the  proper  pore  pressure  boundary  conditions  on  the  upstream  face  of 
the  core. 

The  permeability  of  the  shell  material  in  New  Melones  Dam  is  very 
high  compared  to  that  of  core,  and  therefore  the  excess  pore  pressures 
in  the  shells  of  the  dam  were  assumed  to  be  zero  at  all  stages  of  the 
analysis. 

The  analysis  traced  the  stresses,  the  strains  and  the  pore  pressures 
within  the  dam  throughout  construction  and  for  a period  of  80  years  after 
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the  end  of  construction.  Three  principal  stages  were  considered  in  the 
analysis.  They  are:  the  construction  stage,  the  reservoir  filling  stage, 

and  the  stage  following  reservoir  filling,  when  long-term  seepage  will  be 
established. 

Construction  Stage 

Elements  1 through  16  represent  the  cofferdam.  The  construction  of 
the  cofferdam  was  simulated  in  two  layers  in  this  analysis.  The  con- 
struction of  the  remainder  of  the  dam,  to  an  elevation  of  1135  ft.,  was 
simulated  by  adding  an  additional  six  layers,  consisting  of  elements  17 
through  62. 

Dissipation  of  excess  pore  pressure  during  construction  was  not 
considered  in  the  analysis,  although  the  finite  element  developed  in 
Chapter  6 can  be  used  to  analyze  this  phenomenon.  Dissipation  of  excess 
pore  pressures  during  construction  was  not  considered  in  the  analysis 
because  the  expected  length  of  time  for  construction  of  the  dam  is  only 
three  years. 

The  permeability  of  the  core  is  in  the  range  from  10-7  cm/sec  to 

— 8 

10  cm/sec.  Because  the  core  is  large,  the  amount  of  dissipation  of 
excess  pore  pressure  during  construction  will  be  insignificant.  The  un- 
drained  condition  during  construction  was  simulated  by  using  a zero  time 
step  in  this  portion  of  the  analysis.  Using  a time  step  of  three  or 
even  six  years  would  have  given  virtually  identical  results. 

Reservoir  Filling  Stage 

After  the  dam  has  been  constructed,  the  reservoir  will  be  filled 
to  elevation  990  ft  (MSL) . In  the  analysis,  it  was  assumed  that  the 
water  will  be  raised  to  elevation  990  ft  in  a short  period  of  time,  and 
that  no  dissipation  of  excess  pore  pressure  will  occur  during  this  period. 

The  reservoir  water  pressure  was  applied  on  the  interface  between 
the  core  and  the  shell,  at  the  upstream  side  of  the  dam.  Forces  were 
applied  to  the  submerged  part  of  the  shell  to  simulate  buoyancy.  The 
changes  in  load  at  this  stage  are  shown  in  Fig.  7.5. 
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Long  Term  Seepage  Stage 


Possible  fluctuations  in  the  reservoir  water  level  were  neglected, 
and  it  was  assumed  that  the  reservoir  would  remain  at  elevation  990  ft. 
During  establishment  of  a steady  seepage  condition,  the  deformations 
within  the  dam  will  be  influenced  by  the  dissipation  of  the  excess  pore 
pressures  in  the  core,  and  the  seepage  of  the  water  through  the  core. 

Five  time  steps  were  used  in  this  portion  of  the  analysis.  The 
stresses,  strains,  movements,  and  pore  pressures  were  calculated  at  5, 

15,  30,  50  and  80  years  after  raising  of  the  water  level. 

Settlement 

The  settlement  at  elevation  803  ft  are  shown  in  Fig.  7.6  and 
Fig.  7.7  for  the  two  conditions  considered  in  the  analysis — the  stiff 
core  condition  and  the  soft  core  condition.  Initially,  the  maximum 
settlement,  which  occurs  at  the  center  of  the  core,  is  about  11.5  ft  for 
both  the  stiff  and  the  soft  core.  During  reservoir  filling,  the  core 
deflects  upward  about  1 ft  for  both  the  stiff  and  the  soft  core  condi- 
tion due  to  the  buoyancy  forces  in  the  shell  and  the  submerged  portion 
of  the  core.  During  long  term  seepage,  the  analysis  indicates  approxi- 
mately 1 ft  of  additional  rebound  in  the  stiff  core  case,  and  approxi- 
mately 1 inch  of  rebound  in  the  soft  core  case.  Because  the  initial 
effective  stress  is  greater  in  stiff  core,  the  rebound  capacity  is  higher. 
The  same  effect  can  also  be  noted  in  Fig.  7.8  and  Fig.  7.9  which  show 
the  settlement  along  a nearly  vertical  line  through  the  center  of  the 
core.  It  may  be  noted  that  the  largest  settlement  occurs  near  the  mid- 
height of  the  dam. 

The  fact  that  these  analyses  indicate  some  rebound  during  reservoir 
filling  and  after  is  due  to  the  fact  that  creep  and  secondary  compression 
were  not  considered.  Some  time-dependent  settlements  will  result  from 
these  effects,  and  the  resulting  settlements  would  be  expected  to  affect 
the  rebound  resulting  from  submergence.  The  net  effect  will  probably  be 
a small  amount  of  settlement  rather  than  rebound. 
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FIG.  7.6  SETTLEMENT  ALONG  ELEV.  803  IN  NEW  MELONES  DAM  (STIFF  CORE) 


SETTLEMENT  ALONG  ELEV.  802  IN  NEW  MELONES  DAM  (SOFT  CORE) 


Settlement  (ft) 


FIG.  7.8  SETTLEMENT  ALONG  LINE  B - B’  IN  NEW  MELONES 
DAM  (STIFF  CORE) 
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Horizontal  Movement 


The  horizontal  movements  along  sections  AA'  and  CC'  are  plotted  in 
Figs.  7.10  through  7.13. 

The  calculations  indicate  that  during  construction,  the  core  will 
bulge  at  midheight.  It  may  be  seen  that  the  maximum  calculated  hori- 
zontal movement  on  the  upstream  side  was  about  2 ft  for  both  the  stiff 
and  the  soft  core,  and  the  maximum  calculated  horizontal  movement  on  the 
downstream  side  was  about  1.4  ft  for  the  stiff  core,  and  1.7  ft  for  the 
soft  core.  The  calculations  indicate  that  the  core  will  deflect  down- 
stream about  2 feet  during  reservoir  filling  for  both  the  stiff  core  and 
the  soft  core  cases.  The  results  indicate  that  approximately  the  changes 
in  the  seepage  forces  during  development  of  long  term  seepage  may  cause 
an  additional  three  inches  of  horizontal  movement  in  the  soft  core  case, 
and  approximately  six  inches  in  the  stiff  core  cases.  It  seems  somewhat 
surprising  that  the  calculated  movements  during  consolidation  and  develop- 
ment of  steady  seepage  are  larger  in  the  stiff  core  case  than  in  the  soft 
core  case.  The  reasons  for  this  result  are  believed  to  be:  (1)  the 

movements  are  controlled  more  by  the  shell  stiffness  than  the  core  stiff- 
ness, and  the  shell  stiffness  was  the  same  in  both  cases,  and  (2)  the 
stiff  core  had  lower  pore  pressures  at  the  end  of  construction,  and  the 
load  changes  associated  with  development  of  steady  seepage  were  therefore 
larger  in  the  stiff  core  case. 

Effective  Stresses 

The  calculated  principal  effective  stresses  at  elevation  860  ft 
are  shown  in  Figs.  7.14  through  7.17  for  both  the  stiff  and  the  soft 
core.  In  both  cases  the  major  principal  effective  stress  in  the  shell 
is  higher  than  in  the  core,  because  the  shell  is  stiffer  than  the  core. 
The  major  principal  effective  stress  in  the  stiff  core  is,  as  would  be 
expected,  higher  than  that  in  the  soft  core.  The  results  indicate  that 
during  reservoir  filling,  both  the  major  and  the  minor  stresses  will  be 
reduced  by  1.5-2  tsf  in  the  upstream  side  of  the  shell  due  to  the  effect 
of  buoyancy.  The  calculated  values  of  major  principal  effective  stress 


96 


Elevation  (ft) 


I 


Horizontal  Movement  (ft) 


FIG.  7.11  HORIZONTAL  MOVEMENT  ALONG  SECTION  AA'  IN 
NEW  MELONES  DAM  (SOFT  CORE) 
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FIG.  7.12  HORIZONTAL  MOVEMENT  ALONG  SECTION  CC'  IN 
NEW  MELONES  DAM  (STIFF  CORE) 
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Elevation  (ft) 


FIG.  7.13  HORIZONTAL  MOVEMENTS  ALONG  SECTION  CC’  IN 
NEW  MELONES  DAM  (SOFT  CORE) 
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FIG.  7.14  MAXIMUM  EFFECTIVE  PRINCIPAL  STRESSES  AT 
ELEV.  860  IN  NEW  MELONES  DAM  (STIFF  CORE) 


Minimum  Effective  Principal  Stress,  cr3'(tsf) 


FIG.  7.15  MINIMUM  EFFECTIVE  PRINCIPAL  STRESSES  AT 
ELEV.  860  IN  NEW  MELONES  DAM  (STIFF  CORE) 
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FIG.  7.16  MAXIMUM  EFFECTIVE  PRINCIPAL  STRESSES  AT  ELEV.  860 
IN  NEW  MELONES  DAM  (SOFT  CORE) 


in  the  core  increased  by  about  0.5  tsf,  and  the  calculated  values  of 
minor  principal  effective  stress  increased  about  2.5  tsf  when  the  water 
pressure  was  applied  to  the  upstream  side  of  the  core.  The  increases  in 
the  calculated  stresses  were  slightly  higher  for  the  stiff  core  case  than 
for  soft  core  case. 

During  development  of  steady  seepage,  the  calculated  stresses  on 
the  upstream  side  of  the  core  decreased  because  of  the  seepage  forces 
and  the  calculated  stresses  on  the  downstream  side  of  the  core  increased 
because  of  the  dissipation  of  the  initial  excess  pore  pressures.  The 
calculated  stress  changes  are  somewhat  larger  in  the  stiff  core  case  them 
in  the  soft  core  case. 

Pore  Pressures 

Calculated  pore  pressure  distributions  are  shown  in  Fig.  7.18  and 
Fig.  7.19  for  the  stiff  core  and  the  soft  core  cases  at  different  stages. 
After  construction,  the  calculated  pore  pressure  at  the  base  of  the  soft 
core  was  found  to  be  about  8 tsf,  and  the  corresponding  value  for  the 
stiff  core  was  about  3 tsf.  The  larger  pore  pressure  in  the  soft  core 
occurs  as  a result  of  its  higher  initial  degree  of  saturation  in  the 
soft  core. 

During  reservoir  filling,  the  calculated  pore  pressure  at  the  base 
increased  to  about  10  tsf  in  the  soft  core  case,  and  to  about  8 tsf  in 
the  stiff  core  case,  due  to  the  water  pressure  forces  on  the  upstream 
face  of  the  core. 

The  calculations  indicate  that  both  cores  will  reach  about  the 
same  pore  pressure  distribution  after  80  years.  The  calculated  varia- 
tions of  pore  pressure  with  time  at  points  A and  B at  the  base  of  the 
core  are  shown  in  Figs.  7.20  and  7.21. 

Fig.  7.22  shows  the  calculated  pore  pressure  distribution  in  the 
core  of  the  cofferdam.  It  may  be  seen  that  the  anticipated  pore  pres- 
sures after  construction  are  small.  After  reservoir  filling,  the  calcu- 
lated pore  pressures  in  the  core  increase  due  to  the  weight  of  the  loads 
from  the  water  in  the  reservoir.  The  pore  pressures  at  the  boundary  of 
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FIG.  7.18  PORE  PRESSURE  DISTRIBUTIONS  IN  THE  CORE  OF  NEW  MELONES  DAM  (SOFT  CORE) 


FIG.  7.19  PORE  PRESSURE  DISTRIBUTIONS  IN  THE  CORE  OF  NEW  MELONES  DAM  (STIFF  CORE) 


* 


FIG.  7.20  PORE  PRESSURE  VARIATION  AT  POINT  A IN 
THE  CORE  OF  NEW  MELONES  DAM 
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FIG.  7.21  PORE  PRESSURE  VARIATION  AT  POINT  B 
IN  THE  CORE  OF  NEW  MELONES  DAM 
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FIG.  7.22  PORE  PRESSURE  DISTRIBUTIONS  IN  THE  CORE  OF  THE  COFFERDAM  (STIFF  CORE) 


the  core  and  the  shell  increase,  corresponding  to  the  water  level  at 
elevation  990  ft.  Subsequently,  as  the  water  begins  to  seep  through  the 
core,  the  calculated  pore  pressures  in  the  core  increase  further.  As 
shown  in  the  Fig.  7.20,  after  80  years,  the  results  indicate  that  a 
nearly  steady  state  condition  will  be  reached. 


SUMMARY 

At  the  beginning  of  this  chapter,  three  questions  were  posed  con- 
cerning the  behavior  of  zoned  dams  during  the  time  when  the  excess  pore 
pressures  in  the  core  dissipate  and  steady  seepage  develops.  Based  on 
the  result  discussed  previously,  it  is  possible  to  provide  answers  to 
these  questions  for  conditions  similar  to  these  considered  for  New  Melones 
Dam. 

The  questions,  and  the  answers  indicated  by  the  analyses  of 
New  Melones  Dam,  are: 

(1)  What  is  the  nature  of  the  expected  movements  in  a zoned  dam 
during  consolidation  of  the  core? 

The  maximum  calculated  horizontal  movement  of  the  core  of 
New  Melones  Dam  during  development  of  steady  seepage  was  about  0.6  feet 
downstream.  The  analysis  also  indicates  that  some  upward  movement  may 
occur  during  long  term  seepage  due  to  the  effects  of  buoyancy  as  the  up- 
stream shell  is  submerged.  However,  other  effects  such  as  creep  and 
secondary  compression,  which  were  not  considered  in  the  analysis,  would 
be  expected  to  cause  a net  settlement  rather  than  an  upward  movement 
rebound. 

(2)  How  do  the  stresses  in  the  dam  change  during  consolidation? 

Both  the  maximum  and  minimum  principle  effective  stresses  in  the 
upstream  side  of  the  core  of  New  Melones  Dam  decreased  due  to  the  effects 
of  buoyancy  when  the  upstream  shell  was  submerged  and  long  term  seepage 
developed.  The  maximum  principal  effective  stress  in  the  downstream 
shell  decreased  a small  amount,  and  the  minimum  principal  effective  stress 
increased  a little  during  consolidation  and  development  of  long  term 
seepage  in  the  core. 
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(3)  Are  the  long  term  stresses,  calculated  assuming  slow  con- 
struction and  no  excess  pore  pressure,  the  same  as  these 
calculated  taking  account  of  consolidation? 

For  the  range  of  compaction  condition  considered  likely  for 
New  Melones  Dam,  the  calculated  movements  and  stress  changes  during  con- 
solidation of  the  core  were  quite  small.  The  long-term  stresses  calcu- 
lated by  Quigley,  et  al.  (1976) , using  the  hyperbolic  stress-strain 
relationship  and  assuming  slow  construction  with  no  excess  pore  pres- 
sures at  any  stage,  were  found  to  be  very  close  to  those  calculated  taking 
account  of  consolidation,  as  shown  in  Fig.  7.23. 

Had  a condition  been  analyzed  in  which  the  core  was  much  wetter  and 
softer  than  was  considered  to  be  likely  for  New  Melones  Dam,  it  might  be 
expected  that  the  settlements  and  stress  changes  during  consolidation  of 
the  core  would  be  larger  and  more  important.  It  might  also  be  found 
there  would  be  a greater  difference  between  the  long-term  stresses  calcu- 
lated assuming  slow  construction  with  no  excess  pore  pressures,  and  those 
calculated  taking  account  of  consolidation.  Thus,  although  the  hyper- 
bolic model  and  the  Cam  Clay  Model  produce  very  similar  results  for  the 
case  illustrated  in  Fig.  7.23,  this  conclusion  might  not  apply  to  dams 
with  softer,  wetter  cores. 
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FIG.  7.2  3 COMPARISON  OF  LONG  TERM  STRESS  CALCULATED  USING  HYPERBOLIC  STRESS- 
STRAIN  RELATIONSHIP  WITH  SLOW  CONSTRUCTION  AND  CAM  CLAY  STRFSf- 
STRAIN  RELATIONSHIP  WITH  CONSOLIDATION 


CHAPTER  8 


CONCLUSION 

This  study  was  undertaken  to  develop  a theory  for  analyzing  the 
consolidation  of  partly  saturated  clay  soils,  and  subsequently  to 
develop  a practical  procedure  for  performing  finite  element  analyses  for 
the  behavior  of  dams  during  construction,  during  reservoir  filling,  and 
during  the  development  of  long  term  seepage. 

The  theory  developed  couples  the  effects  of  both  stress  and  flow. 

It  takes  account  of  the  varying  permeability  and  compressibility  of  the 
pore  fluid,  and  the  nonlinear  stress-strain  behavior  of  the  soil. 

In  Chapter  2,  a number  of  elastic-plastic  stress-strain  relation- 
ships were  reviewed  with  respect  to  their  capabilities  for  modelling  the 
stress-strain  behavior  of  compacted  clays.  It  was  found  that  if  the  Cam 
Clay  Model  was  revised  by  introducing  a cohesion  intercept  and  a non- 
linear e-  log  0 ' curve,  it  could  be  used  to  model  the  stress-strain 
m 

behavior  of  compacted  clay  quite  accurately.  This  model  provides  an 
incremental  stress-strain  relationship  which  can  be  used  for  analyses  of 
drained,  undrained,  and  partly  drained  conditions. 

In  Chapter  3,  available  data  were  reviewed  to  determine  how  the 
permeability  of  compacted  clays  varies  with  degree  of  saturation  and  void 
ratio,  and  an  empirical  relationship  was  found  which  expresses  the  depen- 
dence of  permeability  on  these  factors  for  a wide  variety  of  compacted 
clays . 

In  Chapter  4,  a concept  of  homogenized  pore  fluid  was  proposed  to 
simulate  the  air-water  mixture  in  the  pores  of  partly  saturated  clays 
so  that  the  three-phase  partly  saturated  clay  can  be  treated  phenomeno- 
logically as  a two-phase  material.  Compressibility  of  the  homogenized 
pore  fluid  was  then  derived  using  Boyle's  Law  and  Henry's  Law,  and  taking 
account  of  the  effect  of  surface  tension. 
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In  Chapter  5,  Biot's  theory  of  consolidation  was  reviewed,  and  the 
finite  element  formulation  for  the  consolidation  of  saturated  soils  was 
derived.  The  computer  program  based  on  this  formulation  was  verified 
by  comparing  results  derived  from  the  program  with  known  solutions  for 
several  problems. 

In  Chapter  6,  the  theory  of  consolidation  for  partly  saturated  soil 
was  derived  and  a computer  program  was  developed.  To  check  the  utility 
and  effectiveness  of  the  theory  and  the  computer  program  developed,  two 
examples  of  one-dimensional  consolidation  with  compressible  pore  fluid 
and  varying  permeability  were  analyzed. 

In  Chapter  7,  an  analyses  of  New  Melones  Dam  was  performed  using 
the  program  developed.  The  stresses,  strains,  pore  pressures  and  move- 
ments in  the  dam  were  calculated  during  construction  and  at  various 
times  during  a period  of  80  years  following  construction.  To  account 
for  possible  variations  in  the  field  compaction  conditions  in  New  Melones 
Dam,  two  different  conditions  of  water  content  and  compacted  density  of 
the  core  were  considered  in  the  analyses. 

The  analyses  were  performed  for  three  different  stages;  the  con- 
struction stage,  the  reservoir  filling  stage,  and  the  stage  after  < 

reservoir  filling  when  long  term  seepage  had  developed.  The  construction 
and  reservoir  filling  stage  were  assumed  to  occur  under  undrained 
conditions . 

During  the  construction  and  reservoir  filling  stages,  the  pore 
pressures  in  the  soft  core  case  were  found  to  be  significantly  higher 
than  those  in  the  stiff  core  case.  After  a period  of  about  80  years, 
when  a nearly  steady  state  seepage  condition  had  developed,  the  pore 
pressure  distributions  were  found  to  be  about  the  same  for  both  cases. 

More  undrained  deformation  was  found  in  the  soft  core  case  during 
construction  and  reservoir  filling  stages  than  in  the  stiff  core  case, 
and  higher  stresses  were  found  in  the  stiff  core  case  than  in  the  soft 
core  case.  After  reservoir  filling,  the  settlement  and  stress  changes 
during  consolidation  of  the  core  were  small  for  both  the  stiff  core  and 
the  soft  core  cases. 
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Had  a condition  been  analyzed  in  which  the  core  was  much  wetter 
and  softer  than  either  of  the  cases  analyzed,  it  might  be  expected  that 
the  settlements  and  stress  changes  during  consolidation  of  the  core 
would  be  larger  and  more  important  than  either  case  considered  possible 
for  New  Melones  Dam.  It  is  interesting  that  within  a practical  range  of 
conditions  near  optimum  water  content,  the  calculations  indicate  that  the 
movements  and  stress  changes  during  consolidation  can  be  expected  to  be 
quite  small. 

The  studies  conducted  during  this  investigation  have  shown  that  an 
incremental  finite  element  analysis  procedure,  coupled  with  the  non- 
linear stress-strain  behavior,  and  flow  behavior  with  variable  perme- 
ability and  compressibility  of  the  pore  fluid,  is  simple  enough  for 
practical  use,  and  gives  results  which  appear  reasonable  when  compared  to 
field  experience  and  the  results  of  other  analyses.  The  procedure  is 
well  suited  for  practical  use  because  the  required  stress-strain  para- 
meters may  be  determined  from  the  results  of  isotropic  compression  tests 
and  drained  triaxial  tests.  The  other  required  parameters  may  be  readily 
determined  from  compaction  data  and  permeability  test  results. 

Based  on  the  results  of  the  analyses  performed,  it  seems  likely  that 
these  analysis  procedures  may  be  used  to  predict  pore  pressures,  stresses, 
strains  and  movements  in  dams  during  any  stage  after  the  beginning  of 
construction  of  the  dam  and  may  be  useful  for  selecting  desirable  instru- 
ment locations  in  dams,  and  to  help  in  interpreting  the  results  of 
instrumentation  studies 

Perhaps  the  greatest  value  of  these  procedures  is  in  connection 
with  instrumentation  studies.  Piezometers  are  sometimes  used  to  control 
the  rate  of  construction,  and  this  procedure  could  be  even  more  effective 
if  the  rate  of  pore  pressure  development  and  dissipation  could  be  predict- 
ed accurately.  The  analyses  would  also  be  useful  for  interpreting  the 
significance  of  movements  measured  during  and  after  construction. 
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APPENDIX  A 


DERIVATIONS  OF  STRESS-STRAIN  RELATIONSHIP  FOR 
AN  ELASTIC  PLASTIC  MATERIAL 


For  an  elasto-plastic  material  with  hardening  parameter  h,  the 
yield  surface  may  be  expressed  in  the  form 

f({oc’},  h)  = 0 (A. 1) 

For  the  elasto-plastic  model  used  in  this  study,  the  hardening 

parameter  is  a function  of  plastic  strain  only.  {o  '}  is  the  vector  of 

c 

current  effective  stresses.  It  is  related  to  the  changes  in  stress,  dis- 
cussed in  Chapter  5,  as  follows: 

{o  '}  = {o'}  + {o'} 

c u 

where  {o0'}  is  the  vector  of  initial  effective  stresses. 

The  strain  increment  {e}  can  be  thought  of  as  being  composed  of  an 

elastic  component  {e  } and  a plastic  component  {e  },  which  are  related 

E p 

as  follows: 


{e}  = (eE)  + (ep) 

(A. 2) 

where 

{iJ  = [C E]-1{6-> 

E 

(A.  3) 

and  [CE1  is  the  matrix  of  elastic  constants  defined  previously. 

Suppose  that  the  soil  skeleton  has  an  associated  flow  rule 

{y  = x (A-4) 

where  1 is  a positive  parameter,  f is  the  plastic  potential  function. 

For  an  associated  flow  rule,  the  plastic  potential  function  is  identical 
to  the  yield  surface. 

For  an  associated  flow  rule,  substituting  Eqs.  A. 4 and  A. 3 into 
Eq.  A. 2,  may  be  shown  that 
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{e}  = [CE]  1 (or • } + X (A. 5) 

c 

During  plastic  deformation,  the  stress  state  remains  on  the  yield 
surface  so  that 

df  * 3T5iTdW')*!£(wT)d(V  ' 0 <A-6» 

C p 

Using  Eq.  A. 5 and  Eq.  A. 6,  X can  be  expressed  in  terms  of  {e}  as 
follows: 


[cE] {e} 


(aisfr)  lcEl  If  (airr)  (a*7) 


{<5}  = [cE]({e}  - {£  }) 
P 


The  stress-strain  relationship  now  can  be  written  as: 

(A. 8) 

Substituting  Eqs.  A. 4 and  A. 7 into  Eq.  A. 8,  it  may  be  shown  that 

{*}  = [cEp]{e}  (A. 9) 

where  [CEp]  is  the  elasto-plastic  stress-strain  matrix  defined  by 


[CEp]  = [CE]  - 


tcV 


/Jf  \ r^Ei  M /_3h_\  9f 

(rr^f)  3T57T-  8h  („  ,)  9T57T 


(A. 10) 


lcEp)  is  a symmetric  matrix  because  an  associated  flow  rule  has  been 
assumed  in  the  development. 

For  the  Cam  Clay  elasto-plastic  model  described  in  Chapter  2,  the 
yield  surface  is  described  in  two  parts. 
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f. ({o  h.)  - 0 

1 C 1 


(A. 11) 


f_({o  h.)  = 0 

AC  A 


(A. 12) 


Each  of  the  yield  functions  depends  on  a different  hardening  para- 
meter. Within  the  region  defined  by  either  of  these  two  functions, 
equation  A. 10  applies.  However,  at  the  junction  of  two  surfaces,  it  is 
assumed  the  total  plastic  strain  is  the  sum  of  the  contributions  from 
both  of  them  (Koiter,  1953) . 


9f  it 

{y  = A1  dTT^T  + x2  aTcTT 


(A. 13) 


Koiter  has  further  shown  that  the  constants  and  X2  may  be  deter- 
mined from  f = 0,  and  f^  = 0 in  a manner  similar  to  that  for  in  Eq.  A. 7. 
It  follows  that 


*2- 


(STTT-)  lcE'(i) 


t 3fl  \T  ,„E,  3fl  3fl  / 3hl  \ 3fl 

(TOTT)  k 1 5T37T-  STVT 

P 


(A. 14) 


3f, 


and 


Hence 


(wrh) 


! S£2  2T  .J.  5£2  3£2  , Sl>2  xT  S£2 

(JIVT)  ic  1 51571-3^(^77)  3T57T 


(A. 15) 


lc?p]  - IC®]  - 


3£l  \T  -E, 

I 

c 


3«,  . 3h. 


/ JXV  rrEl  1 °ri 

laTvT)  SKTva{e})TO7r 

p 

E 3£2  / 3f2  vT  E 

Ic  1 aTcT/T  (^To^t)  tc^ 

/ 3f2  f f E,  \ *f2  / 3h2  v*  af2 

vaTo^Tj  IC  1 aT^T  " ah2  \3{e  } ) aT^^T 


(A. 16) 
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APPENDIX  B 


ANALYTICAL  SOLUTION  FOR  MATERIAL  WITH  COMPRESSIBLE 
PORE  FLUID  IN  A ONE -DIMENSIONAL  CONSOLIDATION  CASE 


In  the  one-dimensional  consolidation  case,  Biot's  equations  reduce 


9 k 9u 

3z  Y(l)  a* 


_1_  1 9u 

D'  Q 9t 


where  D'  is  the  drained  constrained  modulus  of  the  material,  Q is  the 
compressibility  of  the  pore  fluid  and  p is  pore  pressure. 

If  the  permeability  is  constant,  Eq.  B.l  can  also  be  expressed  as 


— 9zu  _ 9u 
°v  9z2  ” 9t 


where 


=v  Cv  / (1  + q) 


_ kDj 

V Y 


For  a material  with  a compressible  pore  fluid,  a change  in  load  is 
carried  by  both  the  soil  skeleton  and  the  pore  fluid.  Therefore,  there 
is  an  initial  settlement  due  to  the  portion  of  the  load  carried  by  the 
soil  skeleton,  and  there  is  an  initial  excess  pore  pressure  due  to  the 
portion  of  the  load  carried  by  the  pore  fluid. 

The  initial  pore  fluid  is: 


U0‘°z 
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Eq.  B . 2 can  be  solved  for  both  two  way  or  one  way  drainage  boundary 
conditions.  It  is  found  that 


u(z,t)  _ (z  \ 

un  “ FU  ' h2  ) 


(B.7) 


/ z cvt  \ 

where  F^  — , — —)  is  the  same  function  as  given  by  Terzaghi's  solution 

except  that  c in  Terzaghi's  solution  is  replaced  by  another  parameter 
cy,  which  is  defined  by  Equation  B.3.  Thus  Terzaghi's  solution  can  be 
used  to  calculate  the  dissipation  of  excess  pore  pressures  for  materials 
with  compressible  pore  fluids. 

The  total  settlement,  w , at  infinite  time  is: 


w 


a 


(B  .8) 


The  initial  settlement,  wQ,  at  zero  time  is: 


w . = w 
0 a 


1 + n 


(B.9) 


The  degree  of  settlement  can  be  expressed  as: 
w(t)  - w^  , c t 

*0 


woo  - wn  V H2  / 


(B.10) 


where 


/ c t \ 

'Or) 


is  the  s£une  function  as  given  by  Terzaghi  (1924).  There- 
fore, Terzaghi's  solution  can  be  used  to  calculate  the  one-dimensional 
settlement  of  a material  with  compressible  pore  fluid. 
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